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ABSTRACT 


LPI radar is a class of radar systems that possess certain performanee 
characteristics that make them nearly undeteetable by today’s digital intereept receivers. 
This presents a significant tactical problem in the battle space. To detect these types of 
radar, new digital receivers that use sophisticated signal processing techniques are 
required. 

This thesis investigates the use of cyclostationary processing to extract the 
modulation parameters from a variety of eontinuous-wave (CW) low-probability-of- 
intereept (LPI) radar waveforms. The eyclostationary detection techniques deseribed 
exploit the fact that digital signals vary in time with single or multiple periodicities, 
beeause they have spectral correlation, namely, non-zero correlation between certain 
frequeney components, at certain frequency shifts. The use of cyclostationary signal 
processing in a non-cooperative intercept receiver can help identify the partieular emitter 
and can help develop eleetronic attaeks. LPI CW waveforms examined inelude Frank 
eodes, polyphase codes (PI through P4), Frequeney Modulated CW (FMCW), Costas 
frequeneies as well as several frequency-shift-keying/phase-shift-keying (FSK/PSK) 
waveforms. It is shown that for signal-to-noise ratios of OdB and -6 dB, the 
cyclostationary signal processing can extract the modulation parameters necessary in 
order to distinguish among the various types of LPI modulations. 
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EXECUTIVE SUMMARY 


LPI radar is a class of radar systems possessing eertain performance 
eharaeteristies that make today’s digital intercept reeeivers virtually unable to deteet 
them. This presents a significant tactical problem in the battle space. To detect these 
types of radar, the military requires new digital receivers that use sophisticated signal 
processing techniques. 

This thesis investigates the use of cyclostationary processing to extraet the 
modulation parameters from a variety of continuous-wave (CW) low-probability-of- 
intereept (LPI) radar waveforms. The eyclostationary detection techniques described 
exploit the faet that digital signals vary in time with single or multiple periodieity, owing 
to their speetral correlation, namely, non-zero correlation between certain frequency 
components, at eertain frequency shifts. The use of eyelostationary signal processing in a 
non-eooperative intereept reeeiver ean help identify the partieular emitter and can help 
develop eleetronic attacks. LPI CW waveforms examined include Frank codes, polyphase 
eodes (PI, P2, P3 and P4), Frequeney Modulated CW (FMCW), Costas frequeneies as 
well as several frequeney-shift-keying/phase-shift-keying (FSK/PSK) waveforms. This 
thesis shows that for signal-to-noise ratios of OdB and -6 dB, the eyelostationary signal 
processing can extract modulation parameters such as carrier frequency, chip rate, code 
rate and bandwidth, necessary in order to distinguish between the various types of LPI 
radar modulations. 

Two computationally efficient methods of eyelostationary processing were 
implemented: Time Smoothing FFT Aeeumulation and Direet Frequeney Smoothing. It is 
possible to verify that the time smoothing method is more eomputationally efficient than 
the frequeney smoothing for signals with higher eomplexity (polyphase eodes, Frank 
codes, Costas and FSK/PSK). The results from both methods were compared and 
discussed for various LPI modulation types. 

The results due to variations of modulation eharaeteristies are compared and the 
effieiency of both eyelostationary methods for each modulation is measured with relation 
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to the original parameters. This thesis also ineludes eomments on whieh LPI radar signals 
were more suitable for eyelostationary analysis and suggestions for future elassifieation 
systems for these signals, using eombined teehniques. 
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I. INTRODUCTION 


A, LPI RADARS 


Standard surveillance radars always faced the problem of being detected by third- 
party intercept receivers. To mitigate this problem, a low probability of interception has 
become a very common tactical requirement for new radars. The new task of these radars 
is essentially two fold; the radar must meet its specified detection performance with 
minimum radiated power and must measure target characteristics with a modulation that 
is difficult for an intercept receiver to identify. [1] 

Because of their low power, wide BW (BW) and frequency variability, LPI radars 
make detection by modern intercept receivers difficult. The resulting basic definition is 
that LPI radars can detect targets at longer ranges than the modem intercept receivers 
detect the radars. 


This main characteristic of LPI radars may be described by the sensitivity 
advantage (S), defined as the ratio between the signal power needed at the intercept 
receiver to that needed at the LPI radar, which can be expressed mathematically as [2]; 




( 1 . 1 . 1 ) 


where Pm is the power required at the intercept receiver to detect a signal; Prt is the 
power required at the LPI radar receiver to detect a target; <jj is the target’s RCS (Radar 

Cross Section); Gt is the bore sight gain of the LPI radar’s directive transmit antenna; G, 

is the gain of the LPI radar’s transmit antenna side lobe in the direction of the intercept 
receiver; Gi is the gain of the intercept receiver’s antenna; Gr is the radar’s receive 
antenna; Ljr is the loss in the intercept receiver; Lrr is the loss between the radar’s 

antenna and receiver; Rj. is the radar to target range and Rj is the radar to intercept 

receiver range. The sensitivity advantage ((5’) depends on the intercept receiver 
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characteristics and should be a high value, on the order of 50 dB, for a ease where we 
have a simple reeeiver against an LPI radar. 

The sueeess of LPI radars depends on how hard it is for a reeeiver to deteet the 
radars’ emission parameters. The proeessing eapabilities of modern ES (Eleetronie 
Support) equipment are inereasing, leading to more speeifie EPI requirements. On the 
reeeiver side, better results in speetral analysis, for non-eooperative deteetion and 
olassifieation, may be obtained if these radar signals are modeled as eyelostationary. 

All digitally modulated signals are eyelostationary, meaning that their 
probabilistie parameters vary in time with single or multiple periodieities. One property 
that extends from this is that digitally modulated signals have a eertain amount of speetral 
eorrelation. In other words, the signal is eorrelated with frequeney-shifted versions of 
itself (auto-eorrelation) at eertain frequeney shifts. Analyzing LPI radar waveforms using 
eyelostationary modeling is advantageous beeause non-zero eorrelation is exhibited 
between eertain frequeney eomponents when their separation is related to the periodieity 
of interest (e.g., the symbol rate or earner frequeney). The value of that speetral 
separation is referred to as the cycle frequency. 

Two main algorithms stand out as eomputationally effieient tools for 
eyelostationary signal proeessing. The first is the time smoothing Past Pourier Transform 
(PET) Aeeumulation Method and the other is the Direet Prequeney Smoothing Method 
[3]. Both tools are implemented in MATLAB® 6.1 for this thesis. [4] 

B, PRINCIPAL CONTRIBUTIONS 

The objeetive of the researeh deseribed in this thesis was to implement in 
MATEAB® [4], two eomputationally effieient eyelostationary algorithms known as the 
Time Smoothing PET Aeeumulation Method and the Direet Prequeney Smoothing 
Method, defined in [3] and investigate them as an ES reeeiver for proeessing EPI radar 
signals. 

The first step was to generate the EPI signals in a standardized way. The eode 
used was ealled the “EPI Signal Generator,” developed by Pernando Taboada [5], and 
also ineludes eontributions by the author of this thesis. The generated modulations were 
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• BPSK (Binary Phase Shift Keying); 

• FMCW (Frequeney Modulated Continuous Wave); 

• Polyphase Codes (P4, P3, P2 PI and Frank Codes); 

• Costas Codes (Frequeney Hopping - FH); 

• FSK/PSK (Combined Frequency Shift Keying and Phase Shift Keying) with a 
Costas frequency distribution; and 

• FSK/PSK with a target matched frequency distribution. 

Once the signal test matrix was completed, simulations to verify the implementation of 
each algorithm were performed in MATLAB® [4] and the results were compared with 
other receiver signal processing techniques, such as Higher Order Statistics [5], 
Quadrature Mirror Filter Banks [6] and the Wigner Distribution [7]. 

A Graphic User Interface (GUI) was developed in MATLAB® [4] to simplify the 
analysis of the simulation results. The output obtained from the cyclostationary signal 
models were then used to determine the various characteristics of each modulation in 
question. 

Previous work has been done to analyze phase modulation techniques such as 
BPSK (Binary Phase-Shift Keying) and QPSK (Quaternary Phase-Shift Keying) using 
time-smoothing techniques [8, 9]. In this thesis, both frequency and time-smoothing 
techniques are used to analyze various LPI radar modulations and to evaluate the 
measurement of the modulation parameters. 

C. THESIS OUTLINES 

The purpose of this thesis is to document the software implementation of a non- 
cooperative cyclostationary receiver for LPI radar waveforms. The remainder of this 
thesis supports this purpose and is organized as follows. 

Chapter II presents a brief description of Low Probability of Intercept (LPI) 
waveforms and their spectral properties. Two FSK and PSK-combined modulations are 
discussed and analyzed. BPSK, P4, FMCW and Costas Codes are described in depth in 
[5]. 
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Chapter III presents the Cyclostationary signal processing algorithms, a brief 
description of the cyclostationary processor, the MATLAB® [4] tools and the extracted 
parameters description. 

Chapter IV shows the analysis of the different modulation types, their parameters, 
as well as the simulation results. 

Chapter V summarizes the results of this thesis and also makes recommendations 
for future research. 

Appendix A contains the MATLAB® [4] M-files used for implementation of both 
algorithms from [3]. 

Appendix B contains the MATLAB® [4] M-files for the LPI Generator blocks for 
FSK/PSK Costas and Target signals. 

Appendix C contains a table of all LPI radar signals analyzed. 


4 



II. LPI WAVEFORMS DESCRIPTION 


A, BACKGROUND 

LPI radars are especially designed to oppose external third-party receivers that 
attempt to identity the system characteristics and emitter location. Modern electronic 
support (ES) receivers can easily detect the high peak power transmitted by pulsed radars. 
The use of CW (Continuous Wave) modulations and the ability to manage the transmitted 
power limiting emission to the minimum power required to detect typical targets, at the 
required range, make LPI radar signals much less detectable. [5] Besides power 
management, LPI radars modulate their transmissions spreading the energy in frequency 
so that the frequency spectrum of the transmitted signal is wider than required to carry 
the signal’s information. Spreading the signal energy reduces the signal-strength-per- 
information BW. 

LPI waveforms investigated in this thesis include BPSK, LMCW, P4, P3, P2, PI, 
Prank Codes, Costas Codes, PSK/PSK with a Costas frequency distribution, and 
PSK/PSK with a target matched frequency distribution. Refer to Pernando Taboada’s [5] 
thesis work for a detailed description of the other LPI modulations analyzed in this thesis. 
The complete matrix of analyzed signals is shown in Appendix C. This thesis presents the 
analysis of one signal example per modulation type. The analysis of the rest of the signals 
is included in a Technical Report to be published. [10] 

This chapter specifically discusses the two modulation types that combine 
frequency and phase-shift keying (PSK/PSK combined). One modulation type is a 
combination of a Costas frequency-hopping technique and binary phase modulation using 
Barker sequences of different lengths. The second is a frequency-hopping technique that 
uses the characteristic frequency distribution of a desired target, creating a matched PSK, 
which is then modulated with a random-phase keying. 
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B, FSK/PSK COMBINED USING A COSTAS-BASED FREQUENCY¬ 
HOPPING (FH) TECHNIQUE 

This modulation technique is the result of a combination of frequency-shift 
keying based on a Costas frequency-hopping matrix and phase-shift keying using Barker 
sequences of different lengths. A thorough description of the implementation of a Costas 
frequency-hopping technique is in [5]. The purpose of this section is to describe briefly 
the phase encoding applied to a Costas signal, generating the FSK/PSK combined 
waveform. 


In a Costas frequency-hopped signal, the firing order defines what frequencies 
will appear and with what duration. Since we are discussing CW radars, the usual 
terminology does not apply to this case. Instead of a “bursf’ of pulses, we have 
frequencies being continuously emitted during a defined period of time. This period may 
be divided into sub-periods, labeled TV for each frequency. The length of each sub-period 
depends on the sampling interval. During each sub-period, the signal frequency (one of 
the Np frequencies) is modulated by a binary phase sequence according to a Barker 


sequence of length five (-h -h -h - +), seven {+ + + -- + -), eleven (+ + +-+ -- -!--) or 

thirteen For example, the FSK/PSK signal defined by S = 1"^, 

1^ 1^, r, 1^ 2+, 2^ 2+, 2', 2^, 3^ 3^ 3^ 3', 3^, 4^ 4+, 4^ 4', 4+ 5^ 5^ 5^ 5', 5^ 
represents a waveform comprised of Np =5 different frequencies, that are each 


subdivided into five phase slots, labeled TV, according to the Barker sequence of length 
five (+ -h -h - -h). The final waveform may be seen as a binary phase-shifting modulation 
within each frequency hop, resulting in 5 phase slots equally distributed in each 
frequency slot, giving a total of 25 phase slots. 


As illustrated in Figure 1, if we consider Np as the number of frequency hops and 
Np as the number of phase slots of duration TV (Chip Period) in each frequency sub¬ 
period TV, the total number of phase slots in the FSK/PSK waveform is given by: 


N = Np * Np 


( 2 . 1 . 1 ) 
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The block diagram in Figure 1 describes the MATLAB® [4] implementation. The 
user defines which sequence of Costas frequency hops to be used and also how long the 
Barker sequence is (5, 7, 11 or 13). The number of frequency hops is pre-defmed to be 
“seven” and the user may select from two different frequency sequences, varying from 
IkHz to 8kHz. The Costas matrix used in the implementation is the following; 

• Costas Sequence 1—> 4716523 

• Costas Sequence 2—> 2 63875 1. 

Frequencies (kHz) 



Voltage 


+1 

1 

*— Frequency 1 —► 

■*— Frequency 2 —► 


■*— Frequency 5 —► 

Time 

0- 

-1 

+ 

+ 

+ 

- 

+ 


















2 


Tp... 






NpTp= Tp 2NpTp=2Tp (NprOTp NpTp 

b) 

Figure 1 a) FSK/PSK Costas LPI Generator MATLAB® [4] code block 

diagram and, b) general FSK/PSK signal containing Nf frequency 
hops with Np phase slots per frequency. 

The Barker sequence is generated and the frequency-hopping signal is then phase- 
modulated accordingly. For example, if the first Costas sequence is selected, after a phase 
modulation using a Barker sequence of length 5, the final waveform becomes S = 4^, 4^, 
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2^, 3"^, 3^, 3"^, 3', 3^. Figure 2 shows the Power Spectral Density (PSD) plots that reveal 
the spread spectrum characteristic of these signals. The Costas sequence used in the 
following example is always the same and the seven frequency hops are 4, 7, 1, 6, 5, 2, 
and 3kHz . The sampling frequency was 15kHz , satisfying the minimum Nyquist rate 
(fs»2*f) for the largest frequency value. All plots were generated in MATLAB® [4] 
using the routines “fsk_psk_costas.m” and “PAF_FSK_PSK.m”; both listed in Appendix 
B. 



Figure 2 PSD plot for a Costas FH waveform with no phase modulation. 


Figure 3 shows the Costas frequency-hopping waveform time domain, before the 
phase modulation. 

































Time - Seconds 

Figure 3 Time domain plot for a Costas FH waveform with no phase 
modulation. 

Figures 4 to 11 show the PSD, Periodie Ambiguity Function (PAF) and Phase 
Plots for the same frequency sequence as in Figure 2 but now with a phase modulation 
using a Barker-11 sequence. Other signals that may be generated are Costas Frequency- 
Hopping with Barker-5, Barker-7 and Barker-13 phase modulation. These signals are 
going to be analyzed later in Chapter IV. 

The firing orders generated using Costas arrays are designated as optimum in 
reference to the side-lobe behavior of the PAF. As we increase the length of the Barker 
sequence (e.g. from five to eleven), we notice a decrease in the side-lobe level as well. 

The waveforms generated are just examples and various combinations of 
frequency-shift keying and phase-shift keying may be applied to obtain similar results. 
Other examples of these kinds of waveforms are presented in [12]. 
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The PAF plots were performed both for a complete period for all Costas 
frequencies in the sequence and for only one period in one frequency hop. Therefore, the 
Costas PAF characteristics as well as the BPSK PAF characteristics may be compared. 

Figure 4 shows the PSD of the Costas sequence 1 but now with the phase 
modulation using a Barker-11 sequence and with one carrier cycle per phase {cpp). We 
may verify the “flat top” characteristic of the spectrum of LPI signals and also the seven 
“frequency hops.” 



Figure 4 PSD for FSK/PSK Costas FH phase modulated with a Barker-11 
sequence and 1 cpp. 

The phase plot in Figure 5 reveals the Barker-11 sequence phase change. Figure 
6, Figure 7 and Figure 8 show the PAF “thumbtack” [11] characteristic of these types of 
Costas signals. These PAF plots were generated for one period of the whole Costas 
sequence. 
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Figure 5 Phase plot for FSK/PSK Costas FH phase modulated with a Barker- 
11 sequence. 
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Figure 6 PAF contour plot for FSK/PSK Costas FFl phase modulated with a 
Barker-11 sequence (plot of one period for all frequencies in one 
Costas sequence). 
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Normalized z 

Figure 7 PAF delay axis cut for FSK/PSK Costas FFI phase modulated with a 
Barker-11 sequence (plot of one period for all frequencies in one 
Costas sequence). 



Normalized v 

Figure 8 PAF Doppler axis cut for FSK/PSK Costas FH phase modulated with 
a Barker-11 sequence (plot of one period for all frequencies in one 
Costas sequence). 
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Figure 9, Figure 10 and Figure 11 show the modulation for one period of one 
frequency hop of the Costas sequence used in Figure 2. Each frequency hop has similar 
PAF characteristics of a single carrier frequency BPSK signal. The plots of these Figures 
were done for the first carrier frequency of the first Costas sequence. Other modulation 
examples are analyzed later in Chapter IV. 



Figure 9 PAF contour plot for FSK/PSK Costas FH phase modulated with a 
Barker-11 sequence (plot of one period for one frequency in the 
Costas sequence). 
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Figure 10 PAF delay axis cut for FSK/PSK Costas FH phase modulated with a 
Barker-11 sequence (plot of one period for one frequency in the 
Costas sequence). 



Normalized v 


Figure 11 PAF Doppler axis cut for FSK/PSK Costas FH phase modulated with 
a Barker-11 sequence (plot of one period for one frequency in the 
Costas sequence). 
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C. FSK/ PSK COMBINED USING A TARGET-MATCHED FREQUENCY 

HOPPING 

Instead of spreading the energy of the signal equally over a broad BW, this type 
of technique concentrates the signal energy in specific spectral locations of most 
importance for the radar and its typical targets. The signals have a pulse compression 
characteristic, and therefore they can achieve a low probability of intercept. 

The implementation starts with a simulated-target time-radar response. The block 
diagram in Figure 12 describes the signal generation in detail. The target signature data is 
Fourier transformed and the frequency components, their correspondent magnitudes, and 
their initial phases are collected. A random selection process chooses each frequency 
with a probability distribution function defined by the spectral characteristics of the target 
of interest obtained from the Fast Fourier Transform (FFT). That is, the frequencies at the 
spectral peaks of the target (highest magnitudes) are transmitted more often. Each 
“frequency hop,” transmitted during a specific period of time, is also modulated in phase, 
having its initial phase value (from FFT) modified by a pseudo-random spreading-phase 
sequence code of values equally likely to be zero or n radians. [11] The matched 
FSK/PSK radar will then use a correlation receiver with a phase mismatched reference 
signal instead of a perfectly phase matched reference. This allows the radar to generate 
signals that can match a target’s spectral response in both magnitude and phase. [11] 
Only a single target test signal is generated and serves our purpose of testing the 
performance of a digital cyclostationary receiver against these kinds of signals. 
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Figure 12 Block diagram of the MATLAB® [4] implementation of FSK/PSK 
target matched waveform. 


Figure 13 shows the 64 complex points target range radar response plot. Figure 14 
reveals the 64 frequency components that will be selected randomly 256 times. Figure 15 
illustrates the frequency firing order, Figure 16 illustrates the PSD and Figures 17, 18 and 
19 illustrate the PAF properties of these signals. Figure 15 shows the histogram of the 64 
frequency components and shows the number of occurrences of each frequency. Note 
that this is similar to the FFT output or probability distribution shown in Figure 14. The 
following figures show one signal example with 5 carrier periods per phase and 256 
frequency hops. Figure 16 shows the PSD plots and reveals the highly spread-spectrum 
characteristics of this type of modulation. Note the noise-like behavior due to the random 
phase modulation. 
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Figure 13 FSK/PSK target 64 complex points radar range simulated response. 



Frequency 

Figure 14 FSK/PSK target frequency probability distribution of 64 frequency 
components. 
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Frequency (Hz) 

Figure 15 FSK/PSK target 64 frequency components histogram with number of 
occurrences per frequency for 256 frequency hops. 



Figure 16 PSD for FSK/PSK target with 64 frequency components, 256 
frequency hops, random phase modulation and 5 cpp. 
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Figure 17 PAF contour plot for FSK/PSK target with 64 frequency components, 
256 frequency hops, random phase modulation and 5 cpp. 
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Figure 18 PAF delay cut for FSK/PSK target with 64 frequency components, 
256 frequency hops, random phase modulation and 5 cpp. 
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Figure 19 PAF Doppler cut for FSK/PSK target with 64 frequency components, 
256 frequency hops, random phase modulation and 5 cpp. 

An extensive discussion regarding the PAF “thumbtack” characteristic of these types of 
waveforms, as shown in the last three Figures, are presented in Donohoe et al in [11], 

In this chapter we discussed the implementation of two complex LPl Radar 
signals using FSK and PSK techniques combined. A brief theoretical and practical 
tutorial on cyclostationarity processing and its implementation using the FFT 
Accumulation Method and the Direct Frequency Smoothing Method is given in the next 
chapter. 
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III. CYCLOSTATIONARY SIGNAL PROCESSING 
ALGORITHMS AND TUTORIAL 


This chapter briefly explains the eyelostationary proeesses, the time-smoothing 
(FAM) and frequeney-smoothing (DFSM) algorithms and how they were implemented. 
A thorough deseription on eyelostationarity and its properties may be found in [12], [13] 
and [14], 


A. CYCLOSTATIONARY THEORY 

The eyelostationary theory for signal proeessing, as deseribed by William A. 
Gardner [12], involves three main properties: 

• Generation of speetral lines by quadratieally transforming a signal; 

• The statistieal property ealled “seeond-order eyelostationarity,” namely the 
periodie fluetuation of the auto-eorrelation funetion with time; 

• The eorrelation property for signal eomponents in distinet speetral bands. 

The eyelostationary attribute, as it is refleeted in the periodieities of the seeond 
order moments of the signal, ean be interpreted in terms of the generation of speetral lines 
from the signal by putting the signal through a quadratie non-linear transformation. This 
property explains the link between the speetral-line generation property and the statistieal 
property ealled “speetral eorrelation”, eorresponding to the eorrelation that exists 
between the random fluetuations of eomponents of the signal residing in distinet speetral 
bands. The eorrelation integral is very important in theoretieal and praetieal applieations 
and may be defined as 

00 

h{x)= ^ f{u)g{x + u)du (3.1.1) 

-00 

Applying an FFT, it forms a Fourier transform pair given by: 

:s{h{x)]=F{s)G^{s) (3.1.2) 
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If f(x) and g(x) are the same funetion, the integral above is normally ealled the 
autocorrelation funetion and ealled cross-correlation if they differ. The autoeorrelation 
funetion is a quadratie transformation of a signal and may be interpreted as a measure of 
the predietability of the signal at time t +t based on knowledge of the signal at time t. 
[13] 

When eonsidering a time series of length T, the autoeorrelation funetion beeomes 
the time-average autoeorrelation funetion given by 



(3.1.3) 


The non-zero eorrelation (seeond-order periodieity) eharaeteristie of a time series x(t) 
exists, in the time domain, if the equation. 



^0 


(3.1.4) 


where a is the eyele frequeney. R“ (r) is the eyelie auto-eorrelation funetion, also known 
as the “time-frequeney limit autoeorrelation funetion”. The derivation of (3.1.4) from 
(3.1.3) using a non-probabilistie approaeh is developed in [10]. Sinee (3.1.4) is a 
generalization of (3.1.3), when « = 0, the DC eomponent of (3.1.4) yields the time- 
average autoeorrelation funetion of (3.1.3). Therefore, the proeess defined by (3.1.4) is 
able to extraet more information from the signal than the proeess defined by (3.1.3). [13] 

It is well known that the PSD may be obtained from the Fourier Transform of the 
autoeorrelation funetion (3.1.3). [14] 

00 

SAf)=\KO)e-‘^’'’dT (3.1.5) 


In the same manner, it is shown in [3] that the Speetral-Correlation Density (SCD), or 
Cyelie-Speetral Density, may also be obtained from the Fourier Transform of the eyelie 
autoeorrelation funetion (3.1.4) 
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(3.1.6) 


= lim - X J/ + - V; f / - 

2 ) 2 ) 

where a is the cycle frequency and: 

T 

X,{f)=\x{u)e-^^’^^^du (3.1.7) 

-T 

which is the Fourier Transform of the time domain signal x(u). The additional variable 
a leads to a two-dimensional representation S" (/), in the bi-frequency plane ox if, a ) 
plane. [12] 

The speetral eorrelation exhibited by eyelostationary or almost cyclostationary 
processes is completely charaeterized by the eyelie speetra (5*“) or charaeterized 
equivalently by the cyclic autocorrelations (7?“). [12] In praetiee, the cyclie-spectral 

density must be estimated beeause the signals being eonsidered are defined over a finite 
time interval (A f), and therefore the cyelic-spectral density eannot be measured exactly. 
Estimates of the cyelie-speetral density ean be obtained via time or frequency-smoothing 
techniques. In this work we will be able to compare both methods when analyzing LPl 
radar signals. 

An estimate of the SCD can be obtained by the time-smoothed eyelie 
periodogram is given by [10]; 

At 
t-\ - 

j s^juj)du (3.1.8) 

At 

t - 

2 

where 


I R-;(T)e 

—X) 


A (u,f) = —X^ \u,f + -]x* \uj-- 

XTrrr rTf ^ r\ I w r\ 


(3.1.9) 


and At = total observation time of the signal, Tw = short-time EFT window length, and; 
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(3.1.10) 


Xj,^(m,/)= I x{u)e 



2 


is the sliding short-time Fourier Transform, and is a viable solution for computing the 
SCD estimations. Using a graphical explanation, in Figure 20, for some signal x(t) the 
frequency components are evaluated over a small time window (sliding FFT time 
length), along the entire observation time interval At. [12] The spectral components 
generated by each short-time Fourier Transform have a resolution. A/ = .In Figure 


20, L is the overlapping factor between each short-time FFT. In order to avoid aliasing 


and cycle leakage on the estimates, the value of L is defined as F < 




k 


At 





i L \ 

I_I 



Short FFT’s 


J 


J 


Figure 20 Pictorial illustration of the estimation of the time-variant spectral 
periodogram (adapted from [12, 17]). 
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Figure 21 shows that the spectral components of each short-time FFT are 
multiplied, still providing the same resolution capability A/ = jZ , for the cyclic- 

/ -'(T 

spectrum estimates. Note that the dummy variable u has been replaced by the time 
instances At each window {Tw), two components centered about some fo and 

separated by some are multiplied together and the resulting sequence of products is 
then integrated over the total time (At), as shown in (3.1.8). 



Figure 21 Sequence of frequency products for each short-time Fourier 
transforms (adapted from [12, 17]). 


The estimation 5*“ (/)« 5*“ (t, /)^, can be made as reliable and accurate as 

desired for any given t and Af, and for all f by making At sufficiently large, provided that 
equation (3.1.4) exists within the interval At and that a substantial amount of smoothing 
is carried out over At, which leads to the Grenander’s Uncertainty Condition At* A/" » 1 
[12]. This Uncertainty Condition means that the observation time (At) must greatly 
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exceed the time window {Tw), which is used to compute the spectral components. A data 
taper window is also used to minimize the effects of cycle and spectral leakage 
(estimation noise), introduced by frequency component side-lobes. [12] 


If we consider the fact that the cycle frequency estimate is Aa » , it results 

that the estimation of some (^, ) represents a very small area on the bi-frequency plane 

as shown in Figure 22 and since one needs a significant number of estimates to represent 
the cyclic spectrum adequately, it follows that obtaining estimates may become very 
computationally demanding. [12] 



Figure 22 Bi-frequency plane, frequency and cycle frequency resolutions on 
detailed area (adapted from [12, 17]). 
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DISCRETE TIME CYCLOSTATIONARY ALGORITHMS 


The computationally efficient algorithms for implementation of time and 
frequency-smoothing techniques are discussed in [3]. These are the FFT Accumulation 
Method (FAM) and the Direct Frequency-Smoothing Method (DFSM) as described 
below. The temporal and spectral smoothing equivalence is also addressed in [12], The 
computationally efficient algorithms for implementation of time and frequency¬ 
smoothing techniques are extensively discussed in [3]. 


1, The Time-Smoothing FFT Accumulation Method (FAM): 

The time-smoothing FFT Accumulation Method was developed to reduce the 
number of computations required to estimate the cyclic spectrum. [3] This technique 
divides the bi-frequency plane into smaller areas called the channel-pair regions and 
computes the estimates a block at a time using the Fast Fourier Transform. Describing the 
estimated time-smoothed periodogram from Equations (3.1.8) and (3.1.9), in discrete 
terms, yields 


Sl{n,k) 


N 
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N-\ 
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n—0 


— X 
N' 


N' 


f A 

n,k + — 

I 2j 




r 

n,k- 

V 



(3.2.1) 


where 

-jlxkn 

X^,(«,A:)= , (3.2.2) 

«-0 

is the Discrete Fourier Transform of x[n], w[n] is the data taper window (e.g. Hamming 
window) and the discrete equivalents of / and a are k and y respectively. Figure 23 
represents a block diagram [13] used in the implementation of this method in 
MATLAB®. [4] 
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Figure 23 FAM block diagram (adapted from [3, 13]). 

The algorithm consists of three basic stages: computation of the complex 
demodulates (divided into data tapering, sliding N' point Fourier transforming and base 
band frequency-downshift translation sections), then computation of the product 
sequences and smoothing of the product sequences. Making a parallel between the 
variables in equations (3.1.8), (3.1.9) and (3.2.1), we have: 


NAME 

CONTINOUS TIME 

DISCRETE TIME 

SCD 

5 ;,«./)., 


Short EFT Size 

Tw 

N' 

Observation Time 

At 

N 

Time 

t 

n 

Frequency 

f 

k 

Cycle Frequency 

a 

Y 

Grenander’s Uncertainty 

Condition 

M = %»1 



Table 1. Comparison of the estimated time-smoothed Periodogram expressed in 
continuous and discrete time. 
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The parameter N represents the total number of discrete instances within the oservatio 
time and represents the number of pints within the discrete short-time sliding FFT. In 
the FAM algorithm, spectral components of a sequence, x[n], are computed using (3.2.2). 
Two components are multiplied (3.2.1) to provide a sample of a cyclic spectrum estimate 
representing a finite area on the bi-frequency plane called a “channel pair region,” as 
shown in Figure 24. There are channel pair regions in the bi-frequency plane. Note 
the sixteen small channel pair regions corresponding to a value of A = 4 in Figure 24. 

N=4 



Figure 24 Division of bi-frequency plane in channel pair regions (adapted from 
[3, 15]). 

A sequence of samples for any particular area may be obtained by multiplying the same 
two components of a series of consecutive short-time sliding FFT's along the entire 
length of the input sequence. After the channelization performed by an A'-point FFT 
sliding over the data with an overlap of L samples, the outputs of the FFT’s were shifted 
in frequency in order to obtain the complex demodulate sequences (see Figure 23). [3] 
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Instead of computing an average of the product of sequences between the complex 
demodulates, as in (3.1.8), they are Fourier transformed with a P-point (second) FFT. The 

computational efficiency of the algorithm is improved by a factor of L since only 
samples are processed for each point estimate. With the sampling frequency, the cycle 
frequency resolution of the decimated algorithm is defined as (compare to 



), the frequency resolution is 



(compare to Af = 



) and the 


Grenander’s Uncertainty Condition is (compare to At-A/’»1). Figure 25, 

clearly reveals that the estimates toward the top and the bottom (yellow areas) of the 
channel pair region do not satisfy the Uncertainty Condition. 



Figure 25 Cycle frequency and frequency resolutions and the Grenander’s 
Uncertainty Condition (adapted from [3, 13]). 


In order to minimize the variability of these point estimates, we can retain only 


K / 

those cyclic spectrum components that are within of y = ± from the center of the 
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channel pair region. [15] A solution to resolve the entire area of the channel pair region 
without leaving gaps is to apply a data taper window on the frequeney axis to obtain 
better coverage. A Hamming window is applied in this implementation. 


2, Direct Frequency-Smoothing Method: 

Direct frequency-smoothing algorithms first compute the data frequency 
components and then execute spectral-correlation operations direetly on the frequeney 
components. Generally, the Direct Frequency-Smoothing Method (DFSM) is 
computationally superior to indireet algorithms that use related quantities, sueh as the 
Wigner-Ville Distribution. [12] But DFSM is normally less effieient than a time¬ 
smoothing approach. 


Time and frequeney-smoothing equivalenee is discussed in details in [10] and 
derives from the fact that randomly fluctuating statistical spectra (in both time and 

frequeney), sueh as S" (compare to 3.1.8, where 7]^ = and S" 

}if / A/ 


converge in the limit (At ^ GO, A/ O) to the non-random limit spectrum 5*“ (/) if the 

limit autocorrelation (t) exists. Employing discrete frequency averaging and 

following the same analogy as in the time-smoothed cyclic periodogram, the basis for the 
DFSM is the discrete time frequency-smoothed eyclie periodogram represented by 


m 


( 

n,k + — 

A* 

f 

n.k — 

2) 

N 

1 2) 


(3.3.1) 


where 


v-i 


—jlTrkn 


«=0 


(3.3.2) 


is the Discrete Fourier Transform of x\n\, w\n\ is a reetangular window of length N that is 
the total number of points of the FFT related to the total observation time At, y is the 
cycle frequency discrete equivalent, the frequency-smoothing ranges over the interval 

\m\<^^ and AA:« ^ is the frequency resolution diserete equivalent. [12] All 
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considerations made previously in Seetion A (Cyclostationary Theory) are still valid and, 
for M »1, the frequeney and time smoothing methods will yield very similar results. 
The bloek diagram on Figure 26 [13] illustrates the implementation of the “DFSM.m” 
MATLAB® [4] routine. 



Figure 26 DFSM algorithm bloek diagram (adapted from [3, 13]). 

In order to provide full coverage of the bi-frequeney plane at a minimal 
computational expense, Equation (3.3.1) is eomputed along a line of constant cycle 

frequency spacing the point estimates by Ak = ^ method is easier to 

implement and generally is used to validate the time-smoothing “FAM.m” code, but this 
method may become more computationally demanding, mainly in the last bloek in which 
the complex-demodulates product sequences are summed. Considerations on the parallel 
processing of both time and frequency algorithms are diseussed in [13]. Combinations of 
both methods are also advantageous for certain applieations. 

Tthe FAM routine took considerably less time than the DFSM routine. The 
eomputing time for the DFSM was two or three times longer than the FAM for signals 
with a large number of samples. Therefore, for very long signals and large M = AtAf , 
FAM should be used to caleulate the eomplete cyelic speetrum. 


32 



























3. 


GUI Implementation: 


A Graphical User Interface (GUI), as shown in Figure 27, was developed using 
the MATLAB® GUIDE Version 6 [4] in order to facilitate the execution of the 
mathematical model processing. 



Hit <Retum> after 
all the entries bellow: 

1) Enter the data file 
directory (same format); 

2) Enter file name 
(any format); 

3) Enter sampling 
fi'equeney; 

4) Choose the 
frequeney resolution; 

5) Choose the value of 
M; and 

6) Hit the desired 
method button. 


Figure 27 Cyclostationary processing GUI schematic tutorial. 


In this GUI, generated by the MATLAB® [4] routines “cyclo.m”, “FAM.m” and 
“DFSM.m”, the user defines 

• The directory where the data file is to be located; 

• The name of the file (without the extension); 

• The sampling frequency; 

• The desired frequency resolution (); and 

• The Grenander’s Uncertainty Condition M = be chosen to be 

equal to 2 or 4. 

These variables are then used to calculate 
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Desired eyele frequeney resolution: 



N' = 


fs. 


- Number of input ehannels in the ehannelization matrix: 


(rounded to the nearest multiple of 2). For example if 
= 700^ = 875, then N' = 1024; 


Overlap Parameter: the short time FFT’s will be performed and shifted by 
L number of samples for every eomputation (instead of at every sample): 





Number of Columns in the Channelization Matrix: P = 





Total number of points in the input 


proeessed: N = PxL=M ■ 


fs 


data to be 


The results of the simulations and a tutorial on how they should be interpreted 
will be diseussed later in this ehapter. We will also be able to eompare the results of this 
algorithm and the Frequeney-Smoothing approaeh. 


Sinee the number of points to be proeessed N is related to the sampling frequeney, 
the frequeney resolution and the value of M are limited by MATLAB® [4] funetions 
memory eonstraints related to the size of the matrix being handled. The suggested 
optimum value of N, in order to avoid memory errors, may vary between 2,048 and 4,096 
points depending on the signal analyzed and may be given by using the suggested values 
on the table below when entering data in the GUI. 


On the Cyelostationary Proeessing GUI, for an optimum N value of 2048 points 

II 

Sampling Frequeney 

7000 Hz 

15000 Hz 


2 

4 

2 

4 

Frequeney resolution (AA:) 

8 

16 

16 

32 


Table 2. Reeommended variables values for GUI users. 
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The platform for all the simulations was a Pentium IV, 2 GHz CPU clock and 
1GB RAM. The most critical memory constraints appeared while using the “contour.m” 
MATLAB® [4] built-in function, but since this is the plot with the best visualization, we 
decided to keep an average maximum number of points of 2,048 whenever possible. 

Moreover, by increasing the value of N from 2048 to 4096 or more only improved 
the results a little. And since the data files were standardized for comparison with the 
results from the methods described in [5], [6] and [7], the input signals were generated to 
have characteristics that would facilitate measuring for all four processing types 
simultaneously. 

All the files in the Cyclostationary MATLAB® [4] package (“cyclo.m, FAM.m 
and dfsm.m”) can be edited very easily to include new capabilities. The function 
“cyclo.m” triggers the GUI and after entering the data, the user may choose between time 
or frequency-smoothing method by clicking the appropriate button. Additional features 
and options to the GUI (“Cyclo.fig”) may be added by using the MATLAB® GUIDE [4] 
and by editing “cyclo.m” accordingly. 


C. PROCESSING TUTORIAL 


In order to give a general idea of how the results of both processing methods 
should look, this section provides some examples of the signals analyzed. 

1. Test Signals: 


For this example, the input complex signals characteristics are as following; 


Tab 


Carrier Frequency (fc) 

1000 Hz 

Sampling Frequency (fs) 

7000 Hz 


e 3. Test signal characteristics. 


This signal is composed of a single carrier sampled at/j=1000Hz; using the real 
part of the signal we expect that the plot on the bi-frequency plane will have the carrier 
frequency at twice its value, as shown in Figure 28 (DFSM SCD for a Test Signal with N 
= 1024, df= 16 and M= 2). 
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1 

Figure 28 DFSM generated SCD bi-frequency for a test signal (lOOOFIz single 
carrier). 



We see that the frequency content is very simple, except for some estimation noise 
around the highest values of the correlation (red dots). This Spectral-Correlation Density 
was generated using the Frequency-Smoothing Method (DFSM). The test signal in Figure 
29 (FAM SCD for a Test Signal with A = 1024, df= 16 andM= 2) shows that the Time- 
Smoothing Method generates similar results that are a little bit noisier because more 
estimation is involved. 

9 
.8 

7 

.6 
.5 

.4 

.3 
.2 

-3000 -2000 -1000 0 1000 2000 3000 

Cycle Frequency (Hz) 

Figure 29 FAM generated SCD bi-frequency for a test signal (lOOOHz single 
carrier). 
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2 . 


BPSK: 


For the BPSK (Barker Code) example, the signals characteristies are as following: 


Carrier Frequency (fc) 

lOOOHz. 

Sampling Frequency (fs) 

7000Hz. 

Barker Code Length (Np) 

11 bits 

SNR 

Only signal 

Number of Carrier Cycles per Barker Bit (NPBB) 

1 


Table 4. BPSK signal eharacteristies. 

In Figure 30, “k” is the frequeney, “y” is the eyele frequeney, ‘yd” is the carrier. 


‘X” is the sampling frequency, is the code rate = ^TpNp) 


, “Ay>” is the number 


of phase changes, related to the length of the Barker sequence used in the phase 
modulation, and “7p” is the chip period. 
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Figure 30 Pictorial generic illustration of a BPSK signal SCD result. 


37 





















By zooming in on the plot, as shown in Figure 31 and in the examples on Figure 
32 and Figure 33, we can measure the code rate ‘fh \ in the cycle frequency axis. It is 
possible to see a pattern of equally spaced points that will eventually fade away as the 
magnitude decreases (from red to blue). By changing the background color of the contour 
plots in order to facilitate establishing a visual threshold (around 0.4 on the normalized 
values, or “light blue” on the color bar) to eliminate the estimation errors, we can 
approximately determine the BW of the signal being analyzed and this value will then 
divided by the code rate giving the value of (bandwidth = Npxf^^). 



Figure 31 Zoomed in pictorial generic illustration of a BPSK signal SCD result 
with estimation of BW, number of phases and code rate. 
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Figure 32 DFSM generated SCD plot for a BPSK signal (lOOOHz single earner 
and 11 bits Barker-code phase modulation) with estimated BW. 



2320 2340 2360 2380 2400 2420 2440 2460 2480 


Cycle Frequency (Hz) 



Figure 33 Zoomed-in DFSM generated SCD plot for a BPSK signal (estimated 
code rate measurement). 
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The estimation of the BW using this method is somewhat subjeetive and may lead 
to still aeeeptable errors on the order of . From the plots above we may eonsider the 
values: 

• Carrier frequeney ifc) = 1000 Hz (direetly from plot); 

• Code rate 

• BW (direetly from plot) «1000//z ± 90Hz ; 

• Number of Barker Bits (A^) ~ = 11-11« 11 bits. 

This proeedure of zooming-in the plots and measuring the distanees between the 
highest magnitudes (red) points will be used for all the signals analyzed in the next 
ehapter. 


from zoomed-in plots, 


- 1 / 




\TpNp) 


= 90//z 


J 


3. FMCW: 


For the FMCW ease, the signal eharaeteristies are as following: 


Carrier Frequency (/c) 

1000 Hz. 

Sampling Frequency {/$) 

7000 Hz. 

Modulation Band Width ( AF) 

250 Hz 

SNR 

Only signal 

Modulation Period (tm) 

20 ms 


Table 5. FMCW signal eharaeteristies. 

The same behavior showed by BPSK appears here. Figure 34 and Figure 35 show 
that the earrier frequeney appears at twiee its value on the eyele-frequeney axis and that 
the Modulation BW and the modulation period may also be estimated by measuring the 
distanee from the eenter of the SCD result to the edge, on the A: = 0 axis and by measuring 
the distanee (“t/e/to”) between two eonseeutive high eorrelation points in the eyele 
frequeney axis. 
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Figure 34 Pictorial generic illustration of a FMCW signal SCD result. 



Figure 35 Pictorial generic illustration of a zoomed-in plot for a FCMW signal 
SCD. 
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Frequency (Hz) 


This secondary variable called “delta” is related to the modulation period as 
shown in Figure 35 and in the examples of Figure 36 and Figure 37. From Figure 36 
(FAM SCD with N = 2048, M = 4 and AA: = 16) and Figure 37 (Zoomed in plot of 


Figure 36), we get the values: 

• Carrier frequency ifc) = 1000 Flz (directly from plot); 

• BW (AF ± delta ) « 250//z ± 25Hz ; 

• ^in~y2 delta ~ Yl 25~ ~ 2t)ms , where delta=25Hz; 



1700 1800 1900 2000 2100 2200 

Cycle Frequency (Hz) 

Figure 36 FAM generated SCD plot for a FMCW signal (lOOOHz carrier and 
estimated modulation BW of 230 Hz). 
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Figure 37 
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Zoomed-in FAM generated SCD plot for a FMCW signal (“delta ” 
value of 25Hz). 


4. P4: 


The example signal charaeteristics are as following: 


Carrier Frequency (/c) 

1000 Hz 

Sampling Frequency {fs) 

7000 Hz 

Number of Phases {Np) 

16 

SNR 

Only Signal 

Number of Carrier Cycles per Phase (cpp) 

1 


Table 6. P4 signals charaeteristies. 

Figure 38 and Figure 39 show FAM-generated SCD plots. We may see that the 
earner still appears at twiee its value on the cyele frequency axis, where the pattern of 
equally spaced points begins to become visible. Note the white box drawn around the 
data showing the general shape of the SCD plot divided in channel pair regions. In Figure 
39 we can see in detail the estimation of the code rate value {ft). Figure 40 shows the 
estimation of the number of phases, as shown in a pictorial representation on Figure 31. 
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Figure 38 FAM-generated SCD plot for a P4 signal (1 \25Hz carrier and 
estimated BW of 1000//z). 
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Figure 39 Zoomed-in FAM-generated SCD plot for a P4 signal (with estimated 
code rate ifh) of 66Hz). 
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Figure 40 Zoomed-in FAM-generated SCD plot for a P4 signal (with estimation 
of BW and number of phases). 


Therefore 


• Carrier frequency ifc) = 1125Flz (directly from plot); 

• Measured BW 1^5IF = «1000//z ± 66 //z (therefore c/ 7 p=i); 

• Number of Phases Used in the P4 Modulation {Np): 

'^66Hz ~ ~ phases ; 


Np = BW/ =1000//z/ 

/ Jb 


Code Period (7»; Tp = ^PP' = 1' ^^qqq = 0.01 65 ; and 


Code Rate /, - - (% gio) “ 


= 62.5Hz. 
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C. CHAPTER SUMMARY 

This chapter briefly discussed the theory related to the cyclostationarity concept. 
Two computationally efficient cyclostationary processing algorithms (FAM and DFSM) 
were discussed and their implementation in MATLAB® [4] was described. 

The visual feature extraction was explained for some of the modulation types that 
will be analyzed in the next chapter. The contour bi-frequency plots were chosen due to 
their greatest contrast with the figure background and the convenient zooming capability 
of MATLAB® [4] was explored to extraet parameters such as: carrier frequency, code 
rate and code period, BW, number of phases or number of Barker bits present on the 
modulation and modulation period for the FMCW case. In the next section, the signal 
matrix shown in Table 1 is analyzed and the theoretical values are compared with the 
measurements from the SCD plots, according to the tutorial discussed in this chapter. 
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IV. DESCRIPTION OF LPI SPECTRAL PROPERTIES AND 
CYCLOSTATIONARY PROCESSING RESULTS 


In this chapter we will analyze a set of LPI signals, with and without the addition 
of White Gaussian Noise, and extract their main characteristics. All signals were 
generated using the “LPI Generator” code from [5]. A detailed description of all 
waveforms may be found in [5]. Table 8 shows the signal examples analyzed in this 


chapter. 


SIGNAL NAME 

CARRIER OR 
HOPPING 
SEQUENCE 
(kHz) 

BW 

(Hz) 

NPBB/ 

CPP 

NUMBER 
OF BITS / 
PHASES 
AND 
MOD. 
PERIOD 

MOD. TYPE 

B 1 7 7 1 s.mat 

1 

1000 

1 

7 

BPSK 

F 1 7 250 20 s.mat 

1 

250 

X 

20 

FMCW 

PI 1 7 4 1 s.mat 

1 

1000 

1 

16 

PI 

P2 1 7 16 1 s.mat 

1 

1000 

1 

16 

P2 

P3 1 7 16 1 s.mat 

1 

1000 

1 

16 

P3 

P4 1 7 16 1 s.mat 

1 

1000 

1 

16 

P4 

FR 1 7 16 1 s.mat 

1 

1000 

1 

16 

FRANK 

C 1 15 10 s.mat 

4, 7, 1, 6, 5, 2, 3 

- 

10 

1 

COSTAS 

FSK PSK C 1 15 5 1 s.mat 

4, 7, 1, 6, 5, 2, 3 

1000 

1 

5 

FSK/PSK COSTAS 

FSK PSK T 15 128 5 s.mat 

- 

4200 

5 

128 Hops 

FSK/PSK TARGET 

Table 7. List of signal examples analyzec 

1 in this thesis. 


The signal analysis will consist of a brief theoretical description of the signals, an 
illustration of the spectral properties using plots of the PSD, and the resulting plots from 
the Cyclostationary processing algorithms discussed previously. Both time-smoothing 
and frequency-smoothing methods will be used and the method that presents the best 
visualization will be chosen for parameter extraction. The features extracted from the 
signals will be presented in tables comparing the measured results with the given 
parameters. 
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A, TEST SIGNALS 

1, Description 

The Test Signals were chosen to be single carrier signals with no modulation 

• lOOOHz Single Carrier with 7000Hz Sampling Frequency; 

• lOOOHz and 2000Hz Double Carrier with 7000Hz Sampling Frequency. 

These signals give an insight on how the carrier frequency information may be extracted 
from each plot. 

2, Spectral Properties and Results (T_l_7_l_s and T_12_7_2_s) 

Figure 41 shows the PSD of both Test Signals. 




Figure 41 PSD plots for both Test signals; a) lOOOHz single sarrier and b) 1 and 
2000Hz double carrier, no modulation. 
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The single earner ease is diseussed on Chapter III. From Figure 42 and Figure 43 
we ean see that the SCD plots for the real part of the double eanier signal show the inter¬ 
modulation produet at the average of both earner frequeneies (1500Hz) and at 500Hz. 



-6000 -4000 -2000 0 2000 4000 6000 

Cycle Frequency (Hz) 


Figure 42 DFSM generated estimated SCD for a Test signal (lOOOHz and 
2000Hz double earner). 
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1500 2000 2500 3000 3500 4000 4500 
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Figure 43 Zoomed in DFSM generated estimated SCD for a Test signal (1000 
and 2000Hz double eanier). 
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Table 9 compares the measured and original characteristics. 


Feature Extraction 

Characteristic 

Original 

Measured 

Carrier Frequency (/c) 

lOOOHz 

2000Hz 

lOOOHz 

2000Hz 

1500Hz and 500Hz 

(intermodulation 

products) 

Sampling Frequency 

7000Hz. 

given 


Table 8. Comparison between measured and original characteristics for the Test 
signals. 

The Frequency-Smoothing technique generated easier to visualize results, and for such 
simple signals the performance is equivalent to the Time-Smoothing algorithm. A Time- 


Smoothing generated plot is shown in Chapter III for the first Test Signal (lOOOHz). The 
time-smoothing results for the second test signal, with double carriers, are not shown 
since they are very similar to the DFSM plots of Figure 42 and Figure 43. 


B, BPSK 

1, Description: 

Binary Phase Shift Keying codes are not considered LPI radar waveforms in the 
strict sense of the word, but they are here included as a guide or academic tool to 
facilitate the visualization of the resulting outputs from each algorithm. 

This type of signal consists basically of using Barker Sequences of various 
lengths to perform a binary phase modulation of the carrier (0 and n ). The signal is then 
sampled and the output of the generator is saved as two vectors (I and Q) representing the 
in phase and quadrature components of the signal. The modulation is performed 
following the block diagram on Figure 44 [5]. 
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Figure 44 Block diagram for BPSK modulation (from [5]). 


In this work, we are going to analyze a waveform modulated with a Barker 
Sequence of length 7 bits, at various signal-to-noise ratios. All MATLAB® [4] plots were 
generated with the “cyclo.m” generated GUI and using a frequency resolution of 16 Hz 
that influence directly on how precise the measurements of the code rate are made. 

The FAM SCD plots were estimated using A=2048 and M=4, the DFSM SCD plots were 
estimated using A=1024 andM=2. 

2. Spectral Properties and Results (B_l_7_7_l_s) 

The first signal to be analyzed is a 7-Bit Barker Code phase modulated signal, 
with carrier frequency at lOOOHz and with NPBB = 1. No noise is added at this moment 
and the carrier presents a lower magnitude in the PSD plot meaning that this may be 
called a “Suppressed Carrier Dual Side Band” (SC-DSB) Modulation. The first signal 
analyzed has its characteristics defined on Table 10. 


Name 

Carrier 

Sampling 

Number of 

Number of Carrier 


Frequency 

Frequency 

Barker Bits 

Periods per Barker Bit 


(fc) 

{fs) 

{Np) 

(NPBB) 

B 1 7 7 1 s 

1000 Hz 

7000 Hz 

1 

1 


Table 9. B_l_7_7_l_s signal characteristics. 

The following sequence of Figures give an overview of the frequency content (Figure 45) 
and the Estimated SCD (Figure 46, Figure 47, Figure 48, Figure 49 and Figure 50) of the 
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Frequency (Hz) 


signal analyzed. The blue background of the SCD estimation contour plots helps in 
visually filtering out the noise. 



0 1000 2000 :MOO 4000 5000 6000 7000 

Frequency 


Figure 45 PSD for a Barker signal (lOOOHz carrier, 7-bits Barker sequence and 
1 NPBB). 



Cycle Frequency (Hz) 


Figure 46 Estimated FAM SCD contour plot for a BPSK real signal with 
lOOOHz carrier, Barker-7 code and 1 NPBB. 
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Frequency (Hz) 



1000 1500 2000 2500 3000 3500 

Cycle Frequency (Hz) 


Figure 47 Estimated FAM SCD contour plot for a BPSK real signal with 

lOOOHz carrier, Barker-7 code and 1 NPBB, with estimated carrier of 
lOOOHz and estimated BW of lOOOHz. 



2280 2300 2320 2340 2360 2380 2400 2420 2440 
Cycle Frequency (Hz) 


Figure 48 Zoomed-in FAM SCD contour plot for a BPSK signal with lOOOHz 
carrier, Barker-7 code and 1 NPBB, with estimated of 141 Hz. 
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Cycle Frequency (Hz) 


Figure 49 Estimated DFSM SCD eontour plot for a BPSK signal with lOOOHz 
carrier, Barker-7 code and 1 NPBB, with estimated BW of lOOOHz. 
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2250 2300 2350 2400 2450 

Cycle Frequency (Hz) 

Figure 50 Zoomed-in estimated DFSM SCD contour plot for a BPSK signal 
with lOOOHz carrier, Barker-7 code and NPBB = 1, with estimated 
of 142Hz. 
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Frequency (Hz) 


Figures 51 to 55 include the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is included comparing the original and the 
measured characteristics of each signal. 



Figure 51 PSD for a Barker signal (lOOOHz carrier, 7-bits Barker sequence, 1 
NPBB and 0 dB SNR). 



1000 1500 2000 2500 3000 

Cycle Frequency (Hz) 


Figure 52 Estimated FAM SCD contour plot for a BPSK signal with lOOOHz 
carrier, 7-bits Barker code, 1 NPBB and 0 dB SNR, with estimated 
BW of lOOOHz. 
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Figure 53 Zoomed-in estimated FAM SCD for BPSK with lOOOHz carrier, 
Barker-7 code, 1 NPBB, 0 dB SNR, with estimated of 143Hz. 



1000 1500 2000 2500 3000 

Cycle Frequency (Hz) 


Figure 54 Estimated DFSM SCD for BPSK with lOOOFIz carrier, Barker-7 code, 
1 NPBB, 0 dB SNR, with estimated BW of lOOOHz. 
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2250 2300 2350 2400 

Cycle Frequency (Hz) 


Figure 55 Zoomed-in estimated FAM SCD for BPSK with lOOOHz earrier, 
Barker-7 code, 1 NPBB, 0 dB SNR, with estimated of 144Hz. 

The following table shows a comparison between the given and the measured 


parameters. 


Feature Extraction 

Characteristic 

Original 

Measured 

Carrier Frequency (fc) 

1000 Hz 

1000 Hz 

Band Width (BW) 

1000 Hz 

997.5 Hz 

Carrier Periods Per Barker Bit 
(NPBB) 


NPBB = 5 = 1.0025 

Code Rate (/* = 

^100(^ = 142,8 

142.75 Hz (average) 

Number of Barker Bits (Np) 


^"100%42.75" 7-005 


Table 10. Comparison between measured and original characteristics for This BPSK 
signal. 


The BPSK signals analyzed helped as an example in determining how to make 
measurements in more complex signals such as FMCW, Polyphase and frequency 
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hopping signals. Table 11 shows a summary of all measurements for the BPSK 
modulation. Table 12 and Figure 56 show the detection effectiveness of the 


cyclostationary processing for all BPSK cases, comparing with the original values. 


Carrier (Hz) 

Bandwidth (Hz) 

Bits 

Code period (ms) 

SNR 

1000 

990.5 

7 

7.06 

Only signal 

1000 

1004.5 

7 

6.968 

0 

1000 

0 

0 

0 

-6 

1000 

995.5 

11 

11.049 

Only signal 

1000 

1045 

11 

10.526 

0 

1000 

968 

11 

11.363 

-6 

1000 

196 

7 

35.714 

Only signal 

1000 

192.5 

7 

36.363 

0 

1000 

196 

7 

35.714 

-6 

1000 

187 

11 

58.82 

Only signal 

1000 

187 

11 

58.82 

0 

1000 

0 

0 

0 

-6 


Table 11. Summary of all measurements for the BPSK modulation. 


E 

IPSK Detection Efl 

Fectiveness 


Carrier 

Bandwidth 

Code period 

Bits/code 

Only signal 

100% 

98% 

103% 

100% 

OdB 

100% 

99% 

102% 

100% 

(-) 6 dB 

75% 

49% 

51% 

50% 


Table 12. Detection effectiveness for the BPSK modulation. 


BPSK effectiveness on detection 



Parameters 


□ Only signal 

□ OdB 

□ (-)6dB 


Figure 56 Graphic demonstration of detection effectiveness for the BPSK 
modulation. 
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C. FMCW 

1. Description 

The triangular modulation of a Frequency Modulated Continuous Wave (FMCW) 
is very popular as a LPI application for a Continuous Wave signal. The emitter uses a 
continuous 100% duty-cycle waveform and therefore, target range and Doppler 
information can be measured unambiguously while retaining the LPI characteristics. The 
FMCW signal is usually easier to implement than phase code modulation signals as long 
as there is not strict demand on linearity over the modulation BW. [5] The triangular 
modulation consists of two linear frequency modulation sections with positive and 
negative slopes. With this configuration, the range and Doppler frequency of the detected 
target can be extracted unambiguously by taking the sum and the difference of the two 
beat frequencies. These characteristics are showed in Figure 57. 



Figure 57 Linear Frequency Modulated Triangular Waveform and Doppler 
Shifted Signal [5] 
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The frequency of the transmitted signal may calculated by 


^ AF AF 


(4.3.1) 


for 0 < t < tm and zero elsewhere. Here fo stands for RF carrier, AF stands for transmitted 
modulation BW, and tm stands for modulation period. The rate of the frequency change or 
chirp rate F is [5] 


F = 


AF 


(4.3.2) 


The following section shows that we can measure these characteristics of an FMCW 
signal in the bi-frequency plane generated by the cyclostationary algorithms 
implemented. 


2, Spectral Properties and Results (F_l_7_250_20_s) 


This FMCW signal to be analyzed has a modulation period of 20 ms, with carrier 
frequency at lOOOHz and modulation BW of 250 Hz. The signal has its characteristics 
defined on Table 14: 


Name 

Carrier 

Sampling 

Modulation 

Modulation Period 


Frequency 

Frequency 

BW 



ifc) 

ifs) 

(AF) 


F_l_7_250_20_s 

1000 Hz 

7000 Hz 

250 Hz 

20 ms 


Table 13. F_l_7_250_20_s signal characteristics. 


The PSD is shown in Figure 58 and shows the frequency content of the first FMCW 
signal. The following sequence of Figures gives an overview of the frequency content 
(Figure 58), time X frequency plot (Figure 59), and the Estimated SCD (Figure 60, 
Figure 61, Figure 62 and Figure 63, all with X= 2048, frequency resolution of 16 Flz and 
M = 4) of the signal analyzed. 
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Figure 58 PSD for an FMCW signal (lOOOHz carrier, 250Hz modulation BW 
and 20ms modulation period, only signal). 



Figure 59 Time X Frequency plot for an FMCW signal (20ms triangular 
modulation up-ramp period). 
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Figure 60 Estimated FAM SCD eontour plot for an FMCW signal with lOOOHz 
earner and estimated modulation BW of 230Hz. 



2110 2115 2120 2125 2130 2135 2140 

Cycle Frequency (Hz) 



Figure 61 Zoomed-in estimated FAM SCD eontour plot for an FMCW signal 
with an estimated “delta"' value of 25Hz. 
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Frequency (Hz) 



Figure 62 Estimated DFSM SCD contour plot for an FMCW signal with 
lOOOFlz carrier and estimated modulation BW of 235Flz. 



Figure 63 Zoomed-in estimated DFSM SCD contour plot for an FMCW signal 
with an estimated “delta"' value of 21Hz. 
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frequency (Hz) 


Figures 64 to 69 include the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is included comparing the original and the 
measured characteristics of each signal. 



Figure 64 PSD for an FMCW signal (lOOOFlz carrier, 250 Flz modulation BW 
and 20ms modulation period, OdB SNR). 



1800 1900 2000 2100 2200 

cycle frequency (Hz) 

Figure 65 Estimated FAM SCD contour plot for an FMCW signal with lOOOHz 
carrier, 0 dB SNR and estimated modulation BW of 200 Flz. 
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Figure 66 Estimated FAM SCD contour plot for an FMCW signal with an 


estimated “delta ” value of 26Hz. 



Figure 67 PSD for an FMCW signal (lOOOHz carrier, 250 Hz modulation BW 
and 20ms modulation period, -6dB SNR). 
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Frequency (Hz) 



Cycle Frequency (Hz) 

Figure 68 Estimated FAM SCD contour plot for an FMCW signal with lOOOHz 
carrier, -6 dB SNR and estimated modulation BW of 200 Hz. 



2105 2110 2115 2120 2125 2130 2135 


Cycle Frequency (Hz) 

Figure 69 Estimated FAM SCD contour plot for an FMCW signal with an 
estimated “delta ” value of 27Hz. 
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Table 15 shows that the main parameters of the FMCW modulation may be easily 


measured in the bi-frequeney plane, using both time and frequeney-smoothing methods. 


Feature Extraction 

Characteristic 

Original 

Measured 

Carrier Frequency 
(fc) 

1000 Hz 

992.83 Hz (average) 

Modulation BW (AF) 

250 Hz 

218.83 Hz 

^^Delta” 

25 Hz 

24.75 Hz (average) 

Triangular 

Modulation Up-Ramp 
Period (tm) 

C = j 7 ^ = 20ms 

/ 2* delta 

Cl = Yi * 24 75 ~ 0.02020s = 20.20ms 


Table 14. Comparison between measured and original characteristics for an FMCW 


signal. 

Table 16 shows a summary of all measurements for the FMCW modulation. Table 17 and 
graphics on Figure 70 show the detection effectiveness of the cyclostationary processing 


for the FMCW case, comparing with the original values. 


Carrier (Hz) 

Modulation BW (Hz) 

Modulation Period (ms) 

SNR 

992.5 

232.5 

21.74 

Only signal 

993 

208 

19.23 

0 

993 

216 

18.52 

-6 

987.5 

221 

29.42 

Only signal 

987.5 

234 

27.77 

0 

985 

218.72 

0 

-6 

970 

470 

20.83 

Only signal 

970 

460 

17.85 

0 

970 

480 

17.86 

-6 

978 

470 

33.33 

Only signal 

978 

450 

35.71 

0 

970 

480 

0 

-6 


Table 15. Summary of all measurements for the FMCW modulation. 


FMCW Detection Effectiveness 


Carrier 

Modulation Bandwidth 

Modulation Period 

Only signal 

98% 

92% 

106% 

OdB 

98% 

90% 

99% 

(-) 6 dB 

98% 

91% 

45% 


Table 16. Detection effectiveness for the FMCW modulation. 
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FMCW effectiveness on detection 




Parameters 


Figure 70 Graphic demonstration of detection effectiveness for the FMCW 
modulation. 


D. PI 

1, Description 

The first polyphase code to be analyzed is called “PI” and its phase modulation is 
defined by the equation; 


=^VNp-{2j-\mj-wp+ii-m ( 4 . 4 . 1 ) 

Np 

Where i is the number of the samples in a given frequency and j is the index of the 
frequency, the phase of the ith sample of the jth frequency is given by equation (4.4.1), 
and i = \,2,...,Np and j = \,2,...,Np. Therefore, there’s a total of Np^ phases in this 
modulation. Figure 71 reveals the quadratic nature of the phase modulation in PI, given 
by (4.4.1). 
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Figure 71 PI code phase shift. 


2, Spectral Properties and Results (Pl_l_7_16_l_s) 


This PI signal has the following given characteristics: 


Name 

Carrier 

Sampling 

Number of 

Number of Cycles per 


Frequency 

Frequency 

Phases 

Phase {cpp) 


(/c) 

(fs) 

(Np^) 


PI 1 7 16 1 s 

1000 Hz 

7000 Hz 

16 

1 


Table 17. Pl_l_7_16_l_s signal characteristics. 

Based on these properties, we expect a BW of lOOOHz, and yields from (4.4.1) that the 
total number of phases is equal to Np and in this case, the 16 equally spaced regions 


inside the BW are going to be separated by = 



= 100^^ = 62.57/z in the 


cycle frequency axis. The PSD is shown in Figure 72 and shows the frequency content of 
This PI signal. The following sequence of Figures give an overview of the frequency 
content (Figure 72) and the Estimated SCD (Figure 73 and Figure 74, with N = 1024, 
frequency resolution of 16 Hz and M= 2) of the signal analyzed. 
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Frequency (Hz) 



Figure 72 PSD for a PI signal (lOOOHz carrier, 16 phases and 1 cpp, only 
signal). 



Cycle Frequency (Hz) 

Figure 73 Estimated DFSM SCD contour plot for a PI signal with 900Flz carrier 
and estimated BW of lOOOHz. 
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Cycle Frequency (Hz) 

Figure 74 Zoomed-in estimated DFSM SCD eontour plot for a PI signal with an 
estimated code rate ifb) of 62Hz. 

Figures 75 to 79 include the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is included comparing the original and the 
measured characteristics of each signal. 



Figure 75 PSD for a PI signal (lOOOHz carrier, 16 phases and 1 cpp, OdB SNR). 
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Figure 76 Estimated FAM SCD contour plot for a PI signal with 900Hz carrier 
and estimated BW of lOOOHz, with OdB SNR. 
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Cycle Frequency (Hz) 

Figure 77 Zoomed-in estimated FAM SCD contour plot for a PI signal with an 
estimated code rate {ft) of 65Hz, with OdB SNR. 
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78 PSD for a PI signal (lOOOFIz carrier, 16 phases and 
SNR). 


1 cpp, -6dB 
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Figure 79 Estimated FAM SCD contour plot for a PI signal with 850Hz carrier 
and estimated BW of lOOOHz, with -6dB SNR. 
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The FAM measurements eould be performed until an SNR of-6dB, exeept for the 
ff^ value. All results from the DFSM proeessing were similar. 


Feature Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency (/c) 

1000 Hz 

876.67 Hz (average) 

Bandwidth 

(BW) 

1000 Hz 

1000 Hz 

Code Rate (/*) 

62.5 Hz 

63.5 Hz (average) 

Number of 
Phases (Np^) 



Code Period 
itp) 

II 

II 

o 

o 

ON 



Table 18. Comparison between measured and original eharacteristies for a PI signal. 
Both time and frequeney-smoothing proeessing teehniques are very well suited for 
analyzing high BW and short eode period signals. Table 20 shows a summary of all 
measurements for the PI modulation. 


Carrier (Hz) 

Bandwidth (Hz) 

Number of Phases 

Code period (ms) 

SNR 

890 

1000 

4 

16.13 

Only signal 

870 

1000 

4 

15.38 

0 

870 

1000 

0 

0 

-6 

986 

224 

4 

71.43 

Only signal 

986 

210 

4 

71.43 

0 

986 

200 

0 

0 

-6 

986 

1088 

8 

58.82 

Only signal 

975 

990 

8 

66.67 

0 

990 

0 

0 

0 

-6 

1000 

200 

8 

312.5 

Only signal 

1000 

200 

8 

250 

0 

0 

0 

0 

0 

-6 


Table 19. Summary of all measurements for the PI modulation. 


Table 21 and graphies on Figure 80 show the deteetion effeetiveness of the 
eyelostationary proeessing for the PI ease, comparing with the original values for all the 
PI signals analyzed. 
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PI Detection Effectiveness 


Carrier 

Bandwidth 

Code period 

Phases 

Only signal 

97% 

105% 

95% 

100% 

OdB 

96% 

101% 

92% 

100% 

(-) 6 dB 

71% 

50% 

0% 

0% 


Table 20. Detection effectiveness for the PI modulation. 


P1 code effectiveness on detection 



Parameters 


□ Only signal 

□ OdB 

□ (-)6dB 


Figure 80 Graphic demonstration of detection effectiveness for the PI 
modulation. 


E. P2 

1, Description 

The P2 code has almost the same phase increments as the PI code, except that the 
starting phase is different. The P2 code is valid for Np even, and each group of the code 
is symmetric about phase of index 0. These phases may be calculated by the equation 

(I>,j=\^Wp-\)I + (4.5.1) 

Figure 81 reveals the symmetric nature of the P2 code, for an Np^ of 64 phases. 
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Figure 81 P2 symmetric phase relationship between the index in the matrix and 
its phase shift. 


2, Spectral Properties and Results (P2_l_7_16_l_s) 


This P2 signal has the following given properties: 


Name 

Carrier 

Sampling 

Number of 

Number of Cycles per 


Frequency 

Frequency 

Phases 

Phase (cpp) 


(/c) 

ifs) 

(Np^) 


P2 1 7 16 1 s 

1000 Hz 

7000 Hz 

16 

1 


Table 21. P2_l_7_16_l_s signal characteristics. 

From these properties yields a BW of lOOOHz, and from (4.5.1) yields that the total 
number of phases is equal to Np , therefore 16 equally spaced regions inside the BW are 


going to be separated by .2 

/{Np) 



= 62.5Hz in the cycle frequency axis, 


which is related to the code period. The PSD is shown in Figure 82 and shows the 
frequency content of This P2 signal. The following sequence of Figures give an overview 
of the frequency content (Figure 82) and the Estimated SCD (Figure 83 and Figure 84, 
with A^= 1024, frequency resolution of 16 Hz and M= 2) of the signal analyzed. 
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Cycle Frequency (Hz) 

Figure 83 Estimated FAM SCD contour plot for a P2 signal with lOOOHz carrier 
and estimated BW of 950 Hz. 
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Figure 84 Zoomed-in estimated FAM SCD eontour plot for a P2 signal with an 
estimated code rate (ft) of 65 Hz. 

Figures 85 to 87 include the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is included comparing the original and the 
measured characteristics of each signal. 



Figure 85 PSD for a P2 signal (lOOOHz carrier, 16 phases and 1 cpp, OdB SNR). 
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Figure 86 Estimated DSFM SCD contour plot for a P2 signal with lOOOHz 
carrier and estimated BW of 850Hz, with OdB SNR. 
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Figure 87 Zoomed-in estimated DFSM SCD contour plot for a P2 signal with an 
estimated code rate (ft) of 65Hz, with OdB SNR. 
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The results are summarized in Table 23. 


Feature Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency (fc) 

1000 Hz 

1015 Hz (average) 

Bandwidth 

(BW) 

1000 Hz 

900 Hz (average) 

Code Rate (/*) 

62.5 Hz 

65 Hz (average) 

Number of 
Phases {Np^) 


= 13.84 

Code Period {tp) 

II 

II 

o 

o 

ON 



Table 22. Comparison between measured and original eharacteristies for a P2 signal. 


Table 24 shows a summary of all measurements for the P2 modulation. 


Carrier (Hz) 

Bandwidth (Hz) 

Number of Phases 

Code period (ms) 

SNR 

1000 

975 

3.873 

15.38 

Only signal 

1030 

870 

3.66 

15.38 

0 

0 

0 

0 

0 

-6 

1000 

167.5 

3.52 

74.07 

Only signal 

1007 

180 

3.59 

71.43 

0 

1000 

150 

0 

0 

-6 

1000 

962 

8.28 

71.43 

Only signal 

1100 

1100 

4.5 

18.2 

0 

1050 

1200 

0 

0 

-6 

975 

195 

0 

0 

Only signal 

1000 

175 

0 

0 

0 

970 

240 

0 

0 

-6 


Table 23. Summary of all measurements for the P2 modulation. 


Table 25 and graphies on Figure 88 show the deteetion effeetiveness of the 


cyclostationary processing for the P2 case, comparing with the original values. 


P2 Detection Effectiveness 


Carrier 

Bandwidth 

Code period 

Phases 

Only signal 

99% 

94% 

75% 

72% 

OdB 

103% 

94% 

53% 

59% 

(-) 6 dB 

76% 

79% 

0% 

0% 


Table 24. Detection effectiveness for the P2 modulation. 
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P2 code effectiveness on detection 



Parameters 


□ Only signal 

□ OdB 

□ (-)6dB 


Figure 88 Graphic demonstration of detection effectiveness for the P2 
modulation. 


F. P3 

1, Description 

The phase changes in P3 (Figure 89) occur according to the following equation; 

(,- 1)4 

I [(/„+fe)-/,]* = ;rt(,-lf,; (4.6.1) 

0 


where i = 1, 2, .. Np. Substituting k=B/T and tc=l/B, the equation can be written as 

,(P3) _ 7r{i — 1) _ — 


BT 


Np 


(4.6.2) 


where T stands for pulse length, f stands for frequency (/ = f„+kt ),k is a constant, the 
BW is approximately B=kT, tc stands for compressed pulse length {t^ - HB) and the 
provided pulse compression ratio of pc = T ! = BT is provided by the waveform. 
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Figure 89 P3 code phase shift. 


2, Spectral Properties and Results (P3_l_7_16_l_s) 


This P3 signal has the following properties: 


Name 

Carrier 

Sampling 

Number of 

Number of Cycles per 


Frequency 

Frequency 

Phases 

Phase {cpp) 


(/c) 

(fs) 

{Np) 


P3 1 7 16 1 s 

1000 Hz 

7000 Hz 

16 

1 


Table 25. P3_l_7_16_l_s Signal characteristics. 

The expected BW is 1000 Hz, and the total number of phases is equal to 16, therefore 16 
equally spaced regions inside the BW are going to be separated by 


^ ~ cycle frequency axis. The PSD is shown in Figure 

90 and shows the frequency content of This P3 signal. The following sequence of Figures 
give an overview of the frequency content (Figure 90) and the Estimated SCD (Figure 91, 
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Figure 92, Figure 93 and Figure 94, with A^= 1024, frequency resolution of 16 Hz and M 
= 2) of the signal analyzed. 



Figure 90 PSD for a P3 signal (lOOOHz carrier, 16 phases and 1 cpp, only 
signal). 



Cycle Frequency (Hz) 

Figure 91 Estimated FAM SCD contour plot for a P3 signal with 1 lOOHz carrier 
and estimated BW of lOOOHz. 
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Cycle Frequency (Hz) 



Cycle Frequency (Hz) 

Figure 93 Estimated DFSM SCD eontour plot for a P3 signal with 1150Hz 
carrier and estimated BW of lOOOHz. 
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2080 



Cycle Frequency (Hz) 

Figure 94 Zoomed-in estimated FAM SCD contour plot for a P3 signal with an 
estimated ft of 63 Hz. 


Figures 95 to 99 include the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is included comparing the original and the 
measured characteristics of each signal. 



Figure 95 PSD for a P3 signal (lOOOHz carrier, 16 phases, 1 cpp, and OdB 
SNR). 
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500 1000 1500 2000 2500 3000 3500 

Cycle Frequency (Hz) 


Figure 96 Estimated FAM SCD contour plot for a P3 signal with 1050Hz carrier 
and estimated BW of lOOOHz, with OdB SNR. 



2120 2140 2160 2180 2200 

Cycle Frequency (Hz) 

Figure 97 Zoomed-in estimated FAM SCD contour plot for a P3 signal with an 
estimatedyi, of 56 Hz and OdB SNR. 
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igure 98 
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1 

1000 1500 2000 2500 3000 3500 

Cycle Frequency (Hz) 

Estimated DFSM SCD contour plot for a P3 signal with 1050Hz 
carrier and estimated BW of lOOOHz, with OdB SNR. 
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Zoomed-in estimated DFSM SCD contour plot for a P3 signal with an 
estimated ft of 68 Hz and OdB SNR. 
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The given and measured features are summarized in Table 27. 


Feature Extract 

ion 

Characteristic 

Original 

Measured 

Carrier 
Frequency (fc) 

1000 Hz 

1075 Hz 

Bandwidth 

(BW) 

1000 Hz 

1000 Hz 

Code Rate (/*) 

62.5 Hz 

62 Hz (average) 

Numher of 
Phases (Np) 


Np = BW/ ^100^^ = 16 ,3 

Code Period 

II 

II 

o 

o 

ON 

t = ^- = 0.0161295 
p /62 


Table 26. Comparison between measured and original oharaeteristies for a P3 signal. 


Table 28 shows a summary of all measurements for the P3 modulation. 


Carrier (Hz) 

Bandwidth (Hz) 

Numher of Phases 

Code period (ms) 

SNR 

1100 

1000 

16.01 

16 

Only signal 

1050 

1000 

16.03 

16 

0 

0 

0 

0 

0 

-6 

1000 

200 

14.71 

71.43 

Only signal 

1000 

200 

14.71 

71.43 

0 

1000 

220 

0 

0 

-6 

1100 

1000 

71.43 

71.43 

Only signal 

1000 

1000 

71.43 

71.43 

0 

0 

0 

0 

0 

-6 

1000 

200 

7.143 

35.71 

Only signal 

1000 

200 

7.143 

35.71 

0 

0 

0 

0 

0 

-6 


Table 27. Summary of all measurements for the P3 modulation. 


Table 29 and graphics on Figure 100 show the detection effectiveness of the 


cyclostationary processing for the P3 case, comparing with the original values. 


P3 Detection Effectiveness 


Carrier 

Bandwidth 

Code period 

Phases 

Only signal 

105% 

100% 

78% 

79% 

OdB 

101% 

100% 

78% 

79% 

(-) 6 dB 

25% 

28% 

0% 

0% 


Table 28. Detection effectiveness for the P3 modulation. 
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P3 code effectiveness on detection 



Parameters 


□ Only signal 
nOdB 

□ (-) 6dB 


Figure 100 Graphic demonstration of detection effectiveness for the P3 
modulation. 


G. P4 

1, Description 

The P4 code is derived from the same waveform as the P3 code. However, in this 
case, the LO frequency is set equal to f^+kT!2 in the 1-Q detectors. With this 
frequency, the phases of successive samples are defined by the equations; [5] 


( 1 - 1)4 

(l)f^^=2n j [{f^+kt)-f,+kTI2'\dt = 2n j k(t-T/2)dt 


0 


0 


= nk{i - 1)^ - nkT(i - \)t^ 


Np 




(4.7.1) 

(4.7.2) 


Figure 101 shows the relationship between the index in the matrix and its subsequent 
phase for a P4-coded signal with Np= 64 (phases). 
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Figure 101 Phase shift for a P4-eoded signal with A^=64 phases 


2, Spectral Properties and Results (P4_l_7_16_l_s) 

The P4 signal has the following charaeteristies; 


Name 

Carrier 

Sampling 

Number of 

Number of Cycles per 


Frequency 

Frequency 

Phases 

Phase (cpp) 


(/c) 

ifs) 

(Np) 


P4 1 7 16 1 s 

1000 Hz 

7000 Hz 

16 

1 


Table 29. P4_l_7 I6_l_s signal eharacteristies. 


The expeeted BW is 1000 Hz, and the total number of phases is equal to 16. There will be 
16 equally spaeed regions inside the BW separated by ~ ^*^^^6 ~ 

the eycle frequeney axis. The PSD is shown in Figure 102 and shows the frequeney 
eontent of This P4 signal. The following sequenee of Figures give an overview of the 
frequeney eontent (Figure 102) and the Estimated SCD (Figure 103, Figure 104, Figure 
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105, and Figure 106 with N = 1024, frequency resolution of 16 Hz and M = 2) of the 
signal analyzed. 



Figure 102 PSD for a P4 signal (lOOOHz carrier, 16 phases and 1 cpp, only 
signal). 



Cycle Frequency (Hz) 

Figure 103 Estimated FAM SCD contour plot for a P4 signal with 1 lOOHz carrier 
and estimated BW of 1000 Hz. 
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2240 2260 2280 2300 2320 


Cycle Frequency (Hz) 

Figure 104 Zoomed-in estimated FAM SCD eontour plot for a P4 signal with an 
estimated of 66 Hz. 



cycle frequency (Hz) 

Figure 105 Estimated DFSM SCD contour plot for a P4 signal with 1 lOOHz 
carrier and estimated BW of 1000 Hz. 
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Cycle Frequency (Hz) 

Figure 106 Zoomed-in estimated FAM SCD eontour plot for a P4 signal with an 
estimated of 62 Hz. 

Figures 107 to 109 include the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is included comparing the original and the 
measured characteristics of each signal. 



Figure 107 PSD for a P4 signal (lOOOHz carrier, 16 phases, 1 cpp, and OdB 
SNR). 
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Figure 109 Zoomed-in estimated FAM SCD contour plot for a P4 signal with an 
estimated of 67 Hz, with OdB SNR. 
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The given and measured features are summarized in the table below. 


Feature Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency (fc) 

1000 Hz 

1030 Hz (average) 

Bandwidth 

(BW) 

1000 Hz 

1000 Hz 

Code Rate (/*) 

62.5 Hz 

65.5 Hz (average) 

Number of 
Phases (Np) 



Code Period 
(4) 

II 

II 

o 

o 

ON 

C = . = 0.015275 

p /65.5 


Table 30. Comparison between measured and original oharaeteristies for a P4 signal. 


Table 32 shows a summary of all measurements for the P4 modulation. 


Carrier (Hz) 

Bandwidth (Hz) 

Number of Phases 

Code period (ms) 

SNR 

1100 

1000 

15.625 

16 

Only signal 

1100 

1000 

14.93 

16 

0 

1000 

700 

0 

0 

-6 

1000 

200 

14.29 

71.43 

Only signal 

1000 

200 

14.29 

71.43 

0 

960 

120 

0 

0 

-6 

1000 

1000 

71.43 

71.43 

Only signal 

1000 

1000 

71.43 

71.43 

0 

1000 

1100 

0 

0 

-6 

1000 

160 

22.86 

142.86 

Only signal 

1000 

170 

12.59 

74.1 

0 

0 

0 

0 

0 

-6 


Table 31. Summary of all measurements for the P4 modulation. 


Table 33 and graphics on Figure 110 show the detection effectiveness of the 


cyclostationary processing for the P4 case, comparing with the original values. 


P4 Detection Effectiveness 


Carrier 

Bandwidth 

Code period 

Phases 

Only signal 

103% 

95% 

86% 

84% 

OdB 

103% 

96% 

81% 

78% 

(-) 6 dB 

74% 

60% 

0% 

0% 


Table 32. Detection effectiveness for the P4 modulation. 
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P4 code effectiveness on detection 



Parameters 


□ Only signal 

□ 0 dB 

□ (-)6dB 


Figure 110 Graphic demonstration of detection effectiveness for the P4 
modulation. 

H, FRANK 

1, Description 

Frank codes belong to the family of polyphase codes, Chirp codes and Barker 
codes, and have been successfully used in implementing Low probability of Intercept 
(LPI) radar signals. [5] 

■y 

A Frank-coded waveform has a length of Ap and basically consists of a constant 
amplitude signal whose carrier frequency is modulated by the phases of the Frank code 
given by the following equation; 

(4.8.1) 

^ Ap 

where i = number of samples (i = 1,2... Np) and j = number of frequencies (j = 1,2,... 
Ap). Figure 111 clearly displays how the phases change in the Fra nk modulation. 
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2, Spectral Properties and Results (FR_l_7_16_l_s) 


This Frank signal has the following characteristics: 


Name 

Carrier 

Sampling 

Number of 

Number of Cycles per 


Frequency 

Frequency 

Phases 

Phase {cpp) 


(/c) 

(fs) 

{Np^) 


FR 1 7 16 1 s 

1000 Hz 

7000 Hz 

16 

1 


Table 33. FR_l_7_16_l_s signal characteristics. 

The expected BW is 1000 Hz, and the total number of phases is equal to 16. There will be 
sixteen equally spaced regions inside the BW range separated by the distance 


/.= 




= 62.5Hz in the cycle frequency axis. The PSD is shown in 


Figure 112 and shows the frequency content of This Frank signal. The following 
sequence of Figures give an overview of the frequency content (Figure 112) and the 
Estimated SCD (Figure 113 and Figure 114, with N = 1024, frequency resolution of 16 
Hz and M = 2) of the signal analyzed. 
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Frequency (Hz) 



Figure 112 PSD for a Fra nk signal (lOOOHz carrier, 16 phases and 1 cpp, only 
signal). 



500 1000 1500 2000 2500 3000 3500 


Cycle Frequency (Hz) 

Figure 113 Estimated DFSM SCD contour plot for a Fra nk signal with 1150Hz 
carrier and estimated BW of 1000 Hz. 
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Figure 114 Zoomed-in estimated DFSM SCD contour plot for a Fra nk signal with 
an estimatedyi, of 61 Hz. 


Figures 115 to 120 include the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is included comparing the original and the 
measured characteristics of each signal. 



Figure 115 PSD for a Fra nk signal (lOOOHz carrier, 16 phases, 1 cpp, and OdB 
SNR). 
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Figure 116 Estimated FAM SCD contour plot for a Frank signal with lOOOHz 
carrier and estimated BW of 1000 Hz, with OdB SNR. 
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Cycle Frequency (Hz) 

Figure 117 Zoomed-in estimated FAM SCD contour plot for a Fra nk signal with 
an estimated ft of 63 Hz and OdB SNR. 
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Figure 118 
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PSD for a Fra nk signal (lOOOFlz carrier, 16 phases, 1 cpp, and -6dB 
SNR). 
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Figure 119 Estimated FAM SCD contour plot for a Frank signal with 1 lOOHz 
carrier and estimated BW of 1000 Hz, with -6dB SNR. 
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Figure 120 Zoomed-in estimated FAM SCD contour plot for a Fra nk signal with 
an estimated fb of 65 Hz, and -6dB SNR. 




The given and measured features are summarized in Table 35. 


Feature Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency (/c) 

1000 Hz 

1083.3 Hz (average) 

Bandwidth {BW) 

1000 Hz 

1000 Hz 

Code Rate (/*) 

62.5 Hz 

63 Hz (average) 

Number of 
Phases (Np^) 


Np^=BW^^ ^100^^^15,87 

Code Period {tp) 


r = ^, = 0.015875 
p /63 

Table 34. Comparison between measured and origina 

characteristics for a Frank 


signal. 
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Table 36 shows a summary of all measurements for the Frank modulation. 


Carrier (Hz) 

Bandwidth (Hz) 

Number of Phases 

Code period (ms) 

SNR 

1150 

1000 

4.05 

16.4 

Only signal 

1000 

1000 

3.98 

15.87 

0 

1100 

1000 

3.92 

15.38 

-6 

1000 

200 

3.65 

66.67 

Only signal 

1000 

200 

3.65 

66.67 

0 

1000 

210 

0 

0 

-6 

1100 

1000 

7.45 

55.55 

Only signal 

1100 

1000 

8.45 

71.43 

0 

1100 

1000 

0 

0 

-6 

1000 

200 

4.71 

111.11 

Only signal 

1000 

200 

3.92 

76.92 

0 

0 

0 

0 

0 

-6 


Table 35. Summary of all measurements for the Frank modulation. 


The following table and graphies on Figure 121 show the deteetion effectiveness of the 


cyclostationary processing for the Frank case, comparing with the original values. 


Frank Detection Effec 

tiveness 


Carrier 

Bandwidth 

Code period 

Phases 

Only signal 

106% 

100% 

77% 

86% 

OdB 

103% 

100% 

80% 

86% 

(-) 6 dB 

80% 

76% 

24% 

25% 

Table 36. Detection effectiveness for the Frank modu 

ation. 


Frank code effectiveness on detection 



Parameters 


□ Only signal 
nO dB 

□ (-)6dB 


Figure 121 Graphic demonstration of detection effectiveness for the Frank 
modulation. 
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L COSTAS CODES 

1. Description 

We briefly diseussed the Costas modulation on Chapter II. If we don’t phase 
modulate the signals diseussed, we get a pure Costas FH signal. The FH signal eonsists of 
one or more frequeneies being ehosen from a set ...fm) of available frequeneies, 

for transmission at eaeh of a set {t\, t2, ... tn} of eonseeutives time intervals. For 
modeling purposes, it is reasonable to eonsider the situation in whieh m=n, and a 
different one of n equally spaeed frequeneies ...fn) is transmitted during eaeh of 

the n equal duration time intervals {t\, tl, ... tn}. Sueh a signal is represented by a n x n 
permutation matrix A, where the n rows eorrespond to the n frequeneies, the n eolumns 
eorrespond to the n intervals, and the entry ay = 1 means transmission and 0 otherwise. 
[5] Any frequeney sequenee ehosen from a Costas matrix will generate a Costas 
Frequeney Hopping sequenee with the speetral eharaeteristies diseussed in [5]. In order to 
faeilitate the analysis of these signals, the eomplex waveform was analyzed and this will 
result in a less eluttered plot, for this type of signal. 


2. Spectral Properties and Results (C_l_15_10_s) 

The signal eharaeteristies are defined on Table 38. 


Name 

Frequency 

Sequence (kHz) 

Sampling 

Frequency 

ifs) 

Number of Cycles per 

Frequency (cpj) 

C_l_15_10_s 

4716523 

15000 Hz 

10 


Table 37. C_l_15_10_s signal eharaeteristies. 


The expeeted BW is 6000Hz (from lOOOHz to 7000Hz), and the total number of 
frequeneies is equal to 7. The following Figures give an overview of the frequeney 
eontent of the signal (Figure 122) and the Estimated SCD (Figure 123, Figure 124, Figure 
125, Figure 126 and Figure 127, with N= 1024, frequeney resolution of 16 Hz and M = 
2) of the signal analyzed. Note on Figure 123 that the frequeney eomponents are all along 


104 




the / = 0 axis (vertical or frequency axis) while the intermodulation products are located 
symmetrically over the bi-frequency plane, outside the vertical axis. 



Figure 122 PSD for a Costas signal (1, 2, 3, 4, 5, 6 and 7 kHz carriers, 10 cpf, 
only signal). 



Cycle Frequency (Hz) 

Figure 123 Estimated FAM SCD contour plot for a complex Costas signal (1, 2, 
3, 4, 5, 6 and 7000Hz carriers over / = 0 axis, cpj=\0, only signal), 
with intermodulation products. 
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Figures 124 to 127 include the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is included comparing the original and the 
measured characteristics of each signal. 



Figure 124 PSD for a Costas signal (1, 2, 3, 4, 5, 6 and 7kHz carriers, 10 cpf and 
OdB SNR). 


7000 


6000 


5000 


Ivl 

X 


?r'4000 


3000 


2000 


1000 



-6000 -4000 -2000 0 2000 4000 6000 

Cycle Frequency (Hz) 



Figure 125 Estimated FAM SCD contour plot for a complex Costas signal (1, 2, 
3, 4, 5, 6 and 7kHz carriers over 7 = 0 axis, 10 cpf, OdB SNR), with 
intermodulation products. 
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Figure 126 PSD for a Costas signal (1, 2, 3, 4, 5, 6 and 7kHz carriers, 10 cpf and 


SNR of-6dB). 
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Figure 127 Estimated FAM SCD contour plot for a complex Costas signal (1,2, 
3, 4, 5, 6 and 7kHz carriers over 7 = 0 axis, 10 cpf, -6dB SNR), with 
intermodulation products. 
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The order of the frequeneies ean’t be measured in this kind of plot sinee we don’t 


have a '"time x frequency" type of output, but the frequeney eontent is easily determined. 


Feature Extraction 

Characteristic 

Original 

Measured 

Frequency 

Components 

1, 2, 3, 4, 5, 6 and 7kHz 

1, 2, 3, 4, 5, 6 and 7kHz 

BW (BW) 

6000Hz 

6000Hz 


Table 38. Comparison between measured and original eharaeteristies for a Costas 


signal. 


The following table shows a summary of all measurements for the Costas modulation. 


Sequence 

Time in each freq 

Code Period 

SNR 

4716523 

0 

0 

Only signal 

4716523 

0 

0 

0 

4716523 

0 

0 

-6 

4716523 

0 

0 

Only signal 

4716523 

0 

0 

0 

4716523 

0 

0 

-6 

4167532 

0 

0 

Only signal 

4167532 

0 

0 

0 

4167532 

0 

0 

-6 

4167532 

0 

0 

Only signal 

4167532 

0 

0 

0 

4167532 

0 

0 

-6 

Table 39. Summary of all 

. measurements for the Costas modulation. 


The following table and graphies on Figure 128 show the deteetion effeetiveness of the 


eyelostationary proeessing for the Costas ease, eomparing with the original values. 


Costas Detection Effectiveness 


Sequence 

Time in each Frequency 

Code period (ms) 

Only signal 

100% 

0% 

0% 

OdB 

100% 

0% 

0% 

(-) 6 dB 

100% 

0% 

0% 


Table 40. 


Deteetion ef 


feetiveness for the Costas modulation. 
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Costas code effectiveness on detection 



Parameters 


□ Only signal 
nOdB 

□ (-)6dB 


Figure 128 Graphic demonstration of detection effectiveness for the Costas 
modulation. 

J. FSK/ PSK COSTAS 

1, Description 

A detailed description of this type of signal is given in Chapter 11. 


2, Spectral Properties and Results (FSK_PSK_C_l_15_5_l_s) 

This FSK/PSK Costas signal has the following characteristics: 


Name 

Frequency 

Sequence 

(kHz) 

Sampling 

Frequency 

ifs) 

Number of 

Barker 

Bits (Np) 

Number of 

Cycles per 

Phase 

(NPBB) 

FSK PSK C 1 15 5 1 s 

4716523 

15000 Hz 

5 

1 

Table 41. FSK PSK C 1 

15 5 1s signa 

characteristics. 


The expected BW is 7000 Hz (from 500 Hz to 7.5kHz), each frequency will be 
modulated with as BPSK Barker-5 and since the smallest frequency is lOOOHz, the BW 
for each frequency hop will be lOOOHz. Inside each hop, we expect to see a 
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/.= 


BW/ 


1000 / 


'Np ~ /5~ following Figures give an overview of the frequency 

content of the signal (Figure 129), the Barker phase sequence (Figure 130) and the 
Estimated SCD (Figure 131, Figure 132 and Figure 133, with N = 1024, frequency 
resolution of 16 Hz andM= 2) of the signal analyzed. 



Figure 129 PSD for a FSK/PSK Costas signal (1, 2, 3, 4, 5, 6 and 7kHz carriers, 
Barker-5 and 1 NPBB, only signal). 



Figure 130 Barker-5 phase sequence used inside each hop. 
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131 Estimated FAM SCD contour plot for a complex FSK/PSK Costas 
signal (1, 2, 4, 5, 6 and 7kHz measured carriers). 
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Figure 132 Zoomed-in estimated FAM SCD contour plot for a complex 
FSK/PSK Costas signal (4, 5 and 6kHz measured carriers and 
estimated BW of 1000 Hz for each frequency hop). 
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Figure 133 Estimated fb value of 200 Flz for the embedded Barker-5 BPSK 
modulation. 

Figures 134 to 137 inelude the analysis of the same signal with added White Gaussian 
Noise. At the end of this analysis, a table is ineluded comparing the original and the 
measured characteristics of each signal. 
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Figure 134 PSD for a FSK/PSK Costas signal (1, 2, 3, 4, 5, 6 and 7kHz carriers, 
Barker-5 and 1 NPBB, OdB SNR). 
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135 Estimated FAM SCD contour plot for a complex FSK/PSK Costas 


signal (1, 2, 4, 5, 6 and 7kHz measured carriers, OdB SNR). 
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Figure 136 Zoomed-in estimated FAM SCD contour plot for a complex 
FSK/PSK Costas signal (5, 6 and 7kHz measured carriers and 
estimated BW of 1000 Hz for each frequency hop, OdB SNR). 
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Figure 137 Estimated fb value of 200 Hz for the embedded Barker-5 BPSK 
modulation. 


The given and measured features are summarized in the table below. 


Feature Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency (fc) 

1, 2, 3, 4, 5, 6 and 7kHz 

1, 4, 5, 6 and 7kHz 

BW for Each 
Frequency Hop 
(BW) 

1000 Hz 

1000 Hz 

Code Rate (fb) 

200 Hz 

200 Hz 

Number of Bits in 
Barker Code (Np) 


m^-BW/ -1000/ -5 

/fb~ /200 ^ 

Code Period 
Inside Each Hop 

=0,005. 

o 

o 

o 

II 

o 

o 

II 

4 .. 


Table 42. Comparison between measured and original characteristics for a FSK/PSK 
Costas signal. 
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The large BW within each frequency hop may increase the error in the estimation of the 
carrier frequencies. For the example above, it was possible to measure practically all 
characteristics of the signals, with SNR until 0 dB. These signals present a modulation 
complexity level much higher than any modem receiver would be able to analyze. The 
cyclostationary processing algorithms implemented reveal their real processing power 
when facing these types of signals: frequency hopping with large intra-hop BWs and 
small code periods. Both cyclostationary algorithms are very powerful tools to perform 
analysis of LPI signals like the FSK/PSK Costas. Still the biggest constraints reside in 
processing power of the analyzer (computer speed and memory). Table 44 shows a 


summary of all measurements for the Costas modulation. 


Sequence 

Bandwidth 

Np 

Code Period 

SNR 

4 7 1 65 2 

7000 

5 

5 

Only signal 

4 7 1 65 2 

7000 

5 

5 

0 

47 1 

6200 

5 

25 

Only signal 

47 1 

6200 

5 

25 

0 

47 1 

7000 

11 

11 

Only signal 

47 1 

7000 

11 

11 

0 

4 

200 

11 

55 

Only signal 

4 

200 

11 

55 

0 


Table 43. Summary of all measurements for the FSK/PSK Costas modulation. 

Table 45 and graphics on Figure 138 show the detection effectiveness of the 
cyclostationary processing for the FSK/PSK Costas case, comparing with the original 
values. 


FSK/PSK Costas Detection Effectiveness 


Sequence 

Bandwidth (Hz) 

Code period (ms) 

Only signal 

47% 

76% 

100% 

OdB 

47% 

76% 

100% 


Table 44. Detection effectiveness for the FSK/PSK Costas modulation. 
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FSK/PSK Costas effectiveness on detection 
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Sequence Bandwidth (Hz) Code period (ms) 

Parameters 

Figure 138 Graphic demonstration of detection effectiveness for the FSK/PSK 
Costas modulation. 

K. FSK/ PSK TARGET 

1, Description 

Vide Chapter 11 for a detailed description of these types of signals. 


2, Spectral Properties and Results (FSK_PSK_T_15_128_5_s) 

This FSK/PSK Target signal has the following characteristics: 


Name 

Sampling 

Frequency 

ifs) 

Number of 

Frequency Hops 

m 

Number of Cycles 

per Phase (cpp) 

FSK_PSK_T_15_128_5_s 

15000 Hz 

128 

5 


Table 45. FSK_PSK_T_15_128_5_s Signal characteristics. 

The expected BW is around 4200Flz (from 2400Hz to 6600Hz), each frequency will be 
modulated with a random phase sequence of length equal to the number of cycles per 


-^- 1 

□ Only signal 

OOdB 
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phase. The following Figures reveal the frequeney content of the signal (Figure 139), the 
phase plots (Figure 140), the histogram of the frequency hops (Figure 141) and the 
Estimated SCD (Figure 142, with A^= 1024, frequency resolution of 16 Hz and M= 2) of 
the signal analyzed. 



Figure 139 PSD for a FSK/PSK Target signal (4200Hz BW, random phase with 
length 5 and 5 cpp, only signal). 



Figure 140 Random phase sequence of length 5 used inside each hop. 
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Figure 141 Frequency hops histogram after random tiring order generator. 
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Figure 142 Estimated FAM SCD contour plot for a complex FSK/PSK Target 
signal with an estimated BW of 4800Flz. 
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The given and measured features are summarized in the table below. 


Feature Extraction 

Characteristic 

Original 

Measured 

Total BW for Entire Frequency Hopping Sequence 

4200Hz 

4800Hz 


Table 46. Comparison between measured and original oharaeteristies for a FSK/PSK 
Target signal. 


The highest peaks of the speetrum of this signal (as seen in the PSD plot) may also be 
detected in the bi-frequency plane. The firing order of the frequency hopping sequence 
cannot be measured because we don’t have the time information in the output of the 
cyclostationary processing. 

This signal presented the highest degree of complexity of the whole set of signals 
analyzed. The cause of this high complexity may be linked to the presence of two 
pseudo-random factors: 

• The frequency firing order; and 

• The phase modulation for each frequency hop. 

Since the emitter knows the pseudo-random sequence used, it can coherently demodulate 
the signal with little effort but for a non-cooperative receiver this task becomes less 
feasible. Table 48 shows a summary of all measurements for the FSK/PSK Target 
modulation. 


Sequence 

Bandwidth 

Code Period 

SNR 

0 

4800 

0 

Only signal 

0 

3000 

0 

Only signal 

0 

4200 

0 

Only signal 

0 

3000 

0 

Only signal 

Table 47. Summary of all 

. measurements for the FSK/PSK Target modulation. 


Table 49 and graphics on Figure 143 show the detection effectiveness of the 
cyclostationary processing for the FSK/PSK Target case, comparing with the original 
values. The analysis was performed for the only signal case, since the results from the 
added White Gaussian Noise cases were all inconclusive. 
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FSK/PSK Target E 

►election Effectiveness 


Sequence 

Bandwidth 

(Hz) 

Code period 
(ms) 

Only signal 

0% 

95% 

0% 


Table 48. Detection effectiveness for the FSK/PSK Target modulation. 


FSK/PSK Target effectiveness on detection 



□ Only signal 

□ OdB 


Parameters 


Figure 143 Graphic demonstration of detection effectiveness for the FSK/PSK 
Target modulation. 

L. COMPARISON BETWEEN POLYPHASE CODES 

The following figures show a comparison between Frank (Figure 144), PI (Figure 
145), P2 (Figure 146), P3 (Figure 147) and P4 (Figure 148), with the same number of 
phases {l^p or l^p = 16), same carrier frequency (IkHz) and same number of carrier 
cycles per phase {cpp=\'). It is possible to see that there’s some difference between the 
SCD estimated results for the various modulations. The general shape of Frank, P3 and 
P4 are very similar but PI and P2 stand out very easily from the group. The absence of 
time in the plots prevents us to see the phase changes along the code. The number of 
phases may be calculated, as seen on the previous sections, by measuring the carrier 
frequency, the bandwidth and the code period. Both cyclostationary algorithms were very 
effective on determining all parameters for the LPI Radar signals with larger bandwidth 
(cp/>=l) and shorter code periods. 
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Figure 144 Estimated DFSM SCD contour plot for a Fra nk signal with 1150Hz 
carrier, BW of 1000 Hz, Np^=\6 and cpp=\. 
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Figure 145 Estimated DESM SCD contour plot for a PI signal with 900Hz 
carrier, BW of 1000 Hz, Np^=\6 and cpp=\. 
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146 Estimated DFSM SCD contour plot for a P2 signal with 1050Hz 
carrier, BW of 950 Hz, Np^=\6 and cpp=\. 


600 

400 

200 

0 

-200 

-400 

-600 

1000 1500 2000 2500 3000 3500 



Cycle Frequency (Hz) 



Figure 147 Estimated DFSM SCD contour plot for a P3 signal with 1150Hz 
carrier, BW of 1000 Hz, Np=\6 and cpp=\. 
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cycle frequency (Hz) 

Figure 148 Estimated DFSM SCD contour plot for a P4 signal with 1150Hz 
carrier, BW of 1000 Flz, Np=\6 and cpp=\. 


The DFSM plots were chosen just to establish a standard observation but the FAM 
generated plots, as seen on the previous sections, show the same results and take less 
computing time to process. 


M. CHAPTER SUMMARY 

In this last Chapter we analyzed and extracted parameters from a set of FPl Radar 
Signals. The extracted parameters were compared to the original ones in tables at the end 
of every processing. 

The analyzed modulations were 

• BPSK (Binary Phase Shift Keying); 

• FMCW (Frequency Modulated Continuous Wave); 

• Polyphase Codes (P4, P3, P2 PI and Frank Codes); 
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• Costas Codes (Frequency Hopping); 

• FSK/PSK (Combined Frequency Shift Keying and Phase Shift Keying) with a 
Costas frequency distribution; and 

• FSK/PSK (Combined Frequency Shift Keying and Phase Shift Keying) with a 
target matched frequency distribution 

It was possible to see that both time and frequency-smoothing methods of 
cyclostationary processing may be applied to analyze these types of signals yielding 
fairly appropriate results for the majority of the signals. A comparison between the 
polyphase codes is also presented. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


This thesis effort was developed in order to implement two cyclostationary 
processing techniques (Time and Frequency-Smoothing algorithms). The Time 
Smoothing FFT Accumulation Method and the Direct Frequency Smoothing Method 
were implemented in MATLAB® Version 6.1 [4] and tested in simulations using a 
standardized data set. During this work it was possible to verify that the cyclostationary 
processing approach is very suitable for analyzing radar signals with LPI characteristics. 
Depending on the type of modulation, the processing and feature extraction can be more 
or less difficult. The analysis of LPI signals can be improved for any type of signals with 
the combination of the algorithms implemented here and other techniques such as: 

• Parallel Filter Arrays and Higher Order Statistics [5]; 

• Quadrature Mirror Filtering [6]; and 

• Wigner-Ville Distribution [7]. 

The similarities between the DFSM results and the FAM results go until a certain 
level; in the zoomed plots we see that the channel pair regions are a little different in 
shape and size, although they occur in the same values for frequency and cycle 
frequency. This may be the result of the different windows applied in each method 
(Hamming window for FAM and square window for DFSM). The computing time is also 
noticeable larger, two or three times more, for the DFSM routine in comparison to the 
FAM routine. The FAM implementation is recommended for long signals with a large 
number of samples. 

Both methods were very effective for phase modulated signals (BPSK, Frank, PI, 
P2, P3, P4 and FSK/PSK Costas) and for FH signals (Costas) until signal-to-noise ratios 
of OdB and -6dB. The FMCW results showed more variation and errors in the results, 
comparing to the original characteristics of the signals, even on the signal-only cases. The 
FSK/PSK Target signal was also very difficult to analyze but mainly due to its pseudo¬ 
random characteristic and complexity. In the technical report [10], signals with different 
bandwidths and different code rates were included. It is possible to see that the 
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cyclostationary methods are very effeetive for signals with larger bandwidth and larger 
code rates. In order to improve the eomputation time, a parallel proeessing 
implementation is reeommended for future work, using clusters or PC-based parallel 
proeessors for aeademic purposes. 

Starting from the extraeted features, we reeommend a future researeh on the 
automatie elassifieation of these types of signals as a form of maintaining the eontinuity 
of this work. Neural networks have being exhaustively pointed as a viable solution for 
aoeomplishing this task. The use of MATLAB® Version 6.1 and GUIDE Version 2.0 [4] 
was very important to the sueeess of this research and helped the development of a very 
useful tool for future academic applications in Electronic Warfare and Signal Processing 
oriented courses of the Naval Post Graduate Sehool. 
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APPENDIX A. CYCLOSTATIONARY IMPLEMENTATION 
CODES (CYCLO.M, FAM.M AND DFSM.M) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

% CYCLO.M 

% 

% Use: GUI implementation for inputing parameters in Cyelostationary Signal Proeessing Methods: 
DFSM.m and FAM.m 

% 

% Inputs: Path of the file to be analyzed; Name of the file; Sampling Frequeney used in the generation of 
the file to be analyzed, 

% desired frequeney resolution and desired method to be used (DFSM or FAM) 

% 

% Output: Various different plots exploring the bi-frequeney plane representation typieal of eyelostationary 
proeessing 

% 

% Created by Antonio F. Lima, Jr.-lo Ten. Av. Brazilian Air Foree - Mareh/2002 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% These are default eommands and eomments from the MATLAB® GUI guide tool 
% DO NOT EDIT !! 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

funetion varargout = eyelo(varargin) 

% CYCLO Applieation M-fde for eyelo.fig 
% FIG = CYCLO launeh eyelo GUI. 

% CYCLO('eallbaek_name', ...) invoke the named eallbaek. 

% Last Modified by GUIDE v2.0 Ol-May-2002 12:56:05 

if nargin == 0 % LAUNCH GUI 

fig = openfig(mfilename,'reuse'); 

% Use system eolor seheme for figure: 
set(fig,'Color',get(0,'defaultUieontrolBaekgroundColor')); 

% Generate a strueture of handles to pass to eallbaeks, and store it. 
handles = guihandles(fig); 
guidata(fig, handles); 

if nargout > 0 

varargout{l} = fig; 
end 

elseif isehar(varargin{l}) % INVOKE NAMED SUBFUNCTION OR CALLBACK 
try 

if (nargout) 

[varargout{Lnargout}] = feval(varargin{:}); % FEVAL switehyard 

else 

feval(varargin{:}); % FEVAL switehyard 
end 

eateh 
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end 


disp(lasterr); 


end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% These are the GUI functions definitions and global variables 
%- 


function varargout = edit3_Callback(h, eventdata, handles, varargin) 
d = get(h,'String'); 

handles. edits = d; % Save file path in global variable handles.editS (path = handles.editS) 
guidata(h,handles); 


% - 

function varargout = editl_Callback(h, eventdata, handles, varargin) 
name = get(h,'String'); 

handles.editl = name; % Save file name in global variable handles.editl (filename = handles.editl) 

if isstr(name) 
else 

errordlg('Filename: Please enter a valid file name (B_l_7_7_s is default)', 'Bad input', 'modal') 
end 

guidata(h,handles); 


% - 

function varargout = edit2_Callback(h, eventdata, handles, varargin) 
fs = str2double(get(h,'String')); 

handles.edit2 = fs;% Save sampling frequency in global variable handles.edit2 (fs = handles.edit2) 
if isnan(fs) 

errordlg('Sampling Frequency: Please enter a numeric value (Hz)', 'Bad input', 'modal') 
end 

guidata(h,handles); 


% - 

function varargout = popupmenul_Callback(h, eventdata, handles, varargin) 
val = get(h,'Value'); 

switch val % Switches between the values chosen by user on popupmenu 


case 1 

handles.popupmenul = 128;% Save frequency resolution in global variable handles.popupmenuS (df = 
handles .popupmenuS) 

case 2 

handles.popupmenul = 64;% Save frequency resolution in global variable handles.popupmenuS (df = 
handles .popupmenuS) 

case S 

handles.popupmenul = S2;% Save frequency resolution in global variable handles.popupmenuS (df = 
handles .popupmenuS) 
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case 4 

handles.popupmenul = 16;% Save frequency resolution in global variable handles.popupmenu3 (df = 
handles .popupmenu3) 

case 5 

handles.popupmenul = 8;% Save frequency resolution in global variable handles.popupmenu3 (df = 
handles .popupmenu3) 

% 

case 6 

handles.popupmenul = 4;% Save frequency resolution in global variable handles.popupmenu3 (df = 
handles .popupmenu3) 


case 7 

handles.popupmenul = 2;% Save frequency resolution in global variable handles.popupmenu3 (df = 
handles .popupmenu3) 


end 


% More options may be added following the same pattern. 
guidata(h,handles); 


% - 

function varargout = popupmenu2_Callback(h, eventdata, handles, varargin) 

Mval = get(h,'Value'); 

switch Mval % Switches between the values chosen by user on popupmenu 
case 1 

handles.popupmenu2 = 2;% Save reliability condition M in global variable handles.popupmenu2 (M = 
handles .popupmenu2) 

case 2 

handles.popupmenu2 = 4;% Save reliability condition M in global variable handles.popupmenu2 (M = 
handles .popupmenu2) 

case 3 

handles.popupmenu2 = 8;% Save reliability condition M in global variable handles.popupmenu2 (M = 
handles .popupmenu2) 

end 

guidata(h,handles); 


%% - 

function varargout = togglebuttonl_Callback(h, eventdata, handles, varargin) 

dfsm;% Runs frequency smoothing method (DFSM.m), using the values entered - Direct Frequency 
Smoothing Method 


% - 

function varargout = togglebutton2_Callback(h, eventdata, handles, varargin) 

fam;% Runs time smoothing method (FFTAM.m), using the values entered - FFT Accumulation Method 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% FAM.M 

% Cyclostationary Estimator using the time smoothing FFT Accumulation Algorithm. 

% Part of the LPI oriented thesis work with Prof P.E. PACE and Prof Ft. LOOMIS 
% Student: Antonio F. LIMA Jr. - 1st. Lt. Brazilian Air Force 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


% Getting input data from the global variables generated by GUI cyclo s.m and calculating constants 
% present directory = cd; % Save present directory to a variable to make the GUI go back to it in the end 
path = handles.edits; % Saving the path information in a variable 
cd (path);% Loading the fde from the specified path 

fdename = handles.edit 1; % File to be analyzed from GUI generated global variable handles.editl 
fs = handles.edit2; % Sampling frequency from GUI generated global variable handles.edit2 
df = handles.popupmenul; % Frequency Resolution from GUI generated global variable 
handles .popupmenuS 

M = 2; % df/dalpha = M » 1 (reliability condition) 
dalpha = df/M; 

load (filename); %loading file to be analyzed in workspace 

X = 1+j.*Q;% using complex signal (for Real part only, use only I - input data file needs to be formatted as 
% described in (4) 

% Defining parameters 

Np = pow2(nextpow2(fs/df)); % Initially, it is the number of rows in the channelization matrix (X) or 
% the number of input channels defined by fs/df 

L = Np/4; % Overlap factor in order to reduce the number of short time fft's 

% L is the offset between points in the same column at consecutive rows 
% should be less than or equal to Np/4 (Prof. Loomis paper) 

P = pow2(nextpow2(fs/dalpha/L)); % Number of columns formed in the channelization matrix (X) 

N = P*L; % Total number of points in the input data to be processed 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

% Input channelization - this part limits the input total number of points to be analyzed. It also generates a 
Np X P matrix 

% (X) with shifted versions of the input vector in each column. 

if length(x)<N 
x(N) = 0; 

elseif length (x)>N 
x = x(l:N); 
end 


NN = (P-1)*L+Np; 

XX = x; 

xx(NN) = 0; 

XX = xx(:); 

X = zeros(Np,P); 
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for k = 0:P-1 

X(:,k+1) = xx(k*L+l:k*L+Np); % X = (Np X P) Channelization input matrix 

end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Windowing (Hamming window) - a vector of length Np is created with the MATLAB® 'hamming' 
function, then this vector is 

% inserted in a Np X Np matrix diagonal and this result is multiplied by the channelized input matrix (X). 

a = hamming (Np); 

XW = diag(a)*X; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% First FFT - this is the first short time fft (Np points long) applied to each column of matrix X. 

XFl = fft(XW); 

XFl = fftshift(XFl); 

XFl = [XF1(:,P/2+1 :P) XFl(:,l:P/2)]; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%Downconversion - the short sliding fffs results are shifted to base band to obtain decimated complex 
demodulate sequences [2] 

E = zeros(Np,P); 
for k = -Np/2:Np/2-l 
for m = 0:P-1 

E(k+Np/2+l, m+1) = exp(-i*2*pi*k*m*L/Np); 
end 
end 

XD = XF1.*E; 

XD = conj(XD'); % Transposing the matrix and taking the complex conjugate of the signal 
clear ('XFl', 'E', 'XW', 'X', 'x','!', 'Q'); % cleaning up memory space 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Multiplication - the product sequences between each one of the complex demodulates and the complex 
conjugate of the others 

% are formed. This forms the area in the bi-frequency plane as explained in [2] figures 4,7 and 10 

XM = zeros(P,Np^2); 
for k = 1 :Np 
for c = 1 :Np 

XM(:,(k-l)*Np+c) = (XD(:,k).*conj(XD(:,c))); 
end 
end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Second FFT - a P point FFT is applied to XM (in each of its columns) 

XF2 = fft(XM); 

XF2 = fftshift(XF2); 

XF2 = [XF2(:,Np^2/2+l:Np^2) XF2(:,l:Np^2/2)]; 

XF2 = XF2 (P/4:3*P/4,:); 

MM = abs(XF2); 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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%Output matrices range 

alphaO = -fs:fs/N:fs; 
fD = -fs/2:fs/Np:fs/2; 

Sx = zeros(Np+l, 2*N+1); 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Display only the data inside the range of interest - centralizes the bi-frequency plane according to alphaO 
aird fO vectors 

forkl = l:?/2+l 
fork2= l:Np^2 
if rem(k2,Np) == 0 
c = Np/2- 1; 
else 

c = rem(k2,Np) - Np/2 - 1; 
end 

k = ceil(k2/Np)-Np/2-l; 
p = kl-P/4-1; 
alpha = (k-c)/Np-l-(p-l)/N; 
f = (k-l-c)/2/Np; 
ifalpha<-l | alpha>l 
k2=k2-H; 
elseif f<-.5 I f>.5 
k2=k2-H; 
else 

kk= RNp*(f-K.5); 

11 = l-l-N*(alpha -I- 1); 

Sx(rouird(kk), rouird(ll)) = M]Vl(kl,k2); 
eird 
end 
eird 

clear ('alpha', 'XM', 'XF2', 'MM', 'f); % cleaning up memory space 

Sx = Sx./max(max(Sx)); % Normalizes the magnitudes of the values in output matrix (maximum = 1) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

% Plotting Routines 
% saves fde containing results 
% savefile = 'FFTAM result s'; 

% save(savefde,'fdename','N','alphaO', 'fO', 'Sx', 'M', 'df) 

% clear all; 

% clc; 

% 

% load ('FFTAM result s.mat'); 

% figure (df -I- M); 

% mesh 

% view (-55,70); 

% xlabel ('cycle frequency'); ylabel('frequency'); zlabel ('SCD magnitude'); 

% title (['Time Smoothing (FFTAM) generated SCD for ', filename, ', df = ', int2str(df), ' and N = ', 
int2str(N)]); 

% 

% % figure 

% % mesh (alpha0,f0, Sx, 'EdgeColor', 'interp'); 
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% % view (0,90); 

% % xlabel ('cycle frequency'); ylabel('ffequency'); zlabel ('SCD magnitude'); 

% % title (['SCD for ', filename,' N=', int2str(N)]); 

% % % 

% figure 

% subplot(2,2,l:2); 

% plot (fO, Sx); grid; 

% title (['Time Smoothing Processing for ', filename,', df = ', int2str(df),' and N = ', int2str(N)]); 

% xlabel ('frequency (Hz)'); ylabel('SCD magnitude'); 

% subplot(2,2,3:4); 

% plot (alphaO, Sx); grid; 

% xlabel ('cycle frequency (Hz)'); ylabel ('SCD magnitude'); 

% % figure 
%%subplot(2,2,l:2); 

% % contour (alphaO, 10, Sx); grid; 

% % xlabel('cycle frequency (Hz)'); ylabel('frequency (Hz)'); 

% % title (['Time Smoothing (FFTAM) generated SCD for ', filename, ', df = ', int2str(df), ' and N 
int2str(N)]); 

% % zoom(2);colorbar; 

% % subplot(2,2,3:4); 

% % contour (alphaO, 10, Sx); grid; 

% % xlabel('cycle frequency (Hz)'); ylabel('frequency (Hz)'); 

% % zoom(2);coforbar; 

% 

figure 

contour (alphaO, 10, Sx); grid; 

xlabel('cycle frequency (Hz)'); ylabel('frequency (Hz)'); 

title (['Time Smoothing (FFTAM) generated SCD for ', filename, ', df = ', int2str(dl), ' and N 

int2str(N)]); 

colorbar; 

delete ('FFTAM result s.maf); 

% newname = ['FFTAM ' filename]; 

% cd('C:\Documents and Settings\aflima\Tese\Result Figures'); 

% saveas (gcf, newname, 'jpg'); 

cd('C:\Documents and Settings\aflima\Tese\cyclo'); 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% *** Cyclostationary Estimator using the Direct Frequency Smoothing Method algorithm.*** 

% Part of the LPI oriented thesis work with Prof P.E. PACE and Prof Ft. LOOMIS 
% Student: Antonio F. LIMA Jr. - 1st. Lt. Brazilian Air Force 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Getting input data from the global variables generated by GUI cyclo s.m 
% function f = dfsm s (x); 

% present directory = cd; % Save present directory to a variable to make the GUI go back to it in the end 

path = handles.edits; % Saving the path information in a variable 

cd (path);% Loading the fde from the specified path 

filename = handles.editl; % File to be analyzed 

Fs = handles.edit2; % Sampling frequency 

df = handles.popupmenul; % Frequency Resolution 

M = 2; % Smoothing window range 

load (filename); 

X = I;% Using real part of signal 
N = (M*Fs)/df; 

N = pow2 (nextpow2(N)); % windowing record for FFT 


X = fft(x,N); % fft of the truncated (or zero padded) time series 
X = fftshift(X);% shift components of fft 

Xc = conj(X); % precompute the complex conjugate vector 


S = zeros (N,N); 
f = zeros (N,N); 
alfa = zeros (N,N); 
F = Fs/(2*N); 

G = Fs/N; 
m =-M/2+1 :M/2; 


% size of the Spectral Correlation Density matrix 
% size of the frequency matrix; 

% size of the cycle frequency matrix 
% precompute constants - F = Fs/(2*N); 

% precompute constants - G = Fs/N; 

% set frequency smoothing window index 


for k = 1 :N % fix k 

% computes vectors of f and alfa, 

% store frequency and cycle frequency data for given k. 


kl = 1:N; 

f(k,kl) = F*(k+kl-l) - Fs/2; % Computes f values and shift them to center in zero (f = (K+L)/2N) 

[ 1 ] 

alfa(k,kl) = G*(k-kl + N-1) - Fs; % Computes alfa values and shift them to center in zero (alfa = (K- 
L)/N) [1] 


for kl = 1 :N %fix kl = J 

%calculate X(K+m) & conj (X(J+m)) for arguments of X(1:N) only 
B = max(l-k, 1-kl); % Largest min of 1 <= (K+m)| (J+m) <= N 

A = min (N-k, N-kl); % Smallest max of 1 <= (K+m)| (J+m) <= N 
n = m((m<=A) & (m>=B)); %fix the index out of range problem by 

% truncating the window 

if isempty(n) 

S(k,kl) = 0; 
else 

p = k+n; q = kl+n; 

Y = X(p).*Xc(q); 
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S(k,kl) = sum(Y); 
end 
end 
end 

S = abs(S./max(max(S)));% normalize output matrix 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% Saving variables, cleaning up memory and running plots 

savefde = 'DFSMresults'; 

%save(savefde,'r, 'Q', 'alfa', 'f, 'S', 'N', 'Fs', 'fdename', 'M', 'df) 
save(savefde,'alfa', 'f, 'S', 'N','fdename', 'M', 'df) 
clear all; 

load ('DFSM result s.maf); 

% figure (2*(df+M)); 

% mesh (alfa,f, S, 'EdgeColor', 'interp'); 

% view (-55,70); 

% xlabel ('cycle frequency'); ylabel('frequency'); zlabel ('SCD magnitude'); 

% title (['Frequency Smoothing (DFSM) generated SCD for ', filename, ' df = ', int2str(df), ' N = ', 
int2str(N)]); 

% % 

% figure (2*(df+M) -I- 1); 

% mesh (alfa,f, S, 'EdgeColor', 'interp'); 

% view (0,90); 

% xlabel ('cycle frequency'); ylabel('frequency'); zlabel ('SCD magnitude'); 

% title (['SCD for', filename,' df =', int2str(df),' N=', int2str(N)]); 

% 

% figure(2*(df+M) -I- 2) 

%% subplot (2,2,1:2) 

% psd(I, N, Fs); 

% title(['I - real part for nfft = ', int2str(N),' and for default nfft = 256']); 

% subplot (2,2,3:4) 

%psd(I, [],Fs); 

% subplot (2,2,3) 

% psd(I-t-j.*Q, N/2, Fs); 

% title('I-tj.Q - complex'); 

% subplot (2,2,4) 

% psd(I-j.*Q, N/2, Fs); 

% title('I-j.Q - complex'); 

% % 

% figure 

% subplot(2,2,l:2); 

% plot (f, S); grid; 

% title (['Frequency Smoothing Processing for ', filename,', df = ', int2str(df),' and N = ', int2str(N)]); 

% xlabel ('frequency (Flz)'); ylabel('SCD magnitude'); 

% subplot(2,2,3:4); 

% plot (alfa, S); grid; 

% xlabel ('cycle frequency (Flz)'); ylabel ('SCD magnitude'); 

% % % 

% figure (2*(df+M) -I- 5) 

% subplot(2,2,l:2); 

% contour (alfa, f, S); grid; 

% xlabel('cycle frequency (Flz)'); ylabel('frequency (Hz)'); 
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% title (['Frequency Smoothing (DFSM) generated SCD for filename, df = int2str(df), ' and N 
int2str(N)]); 

% zoom(2);colorbar; 

% subplot(2,2,3:4); 

% contour (alfa, f, S); grid; 

% xlabel('cycle frequency (Hz)'); ylabel('frequency (Hz)'); 

% zoom(2);colorbar; 
figure 

contour (alfa, f, S); grid; 

xlabel('cycle frequency (Hz)'); ylabel('frequency (Hz)'); 

title (['Frequency Smoothing (DFSM) generated SCD for ', filename, ', df = ', int2str(df), ' and N 

int2sb(N)]); 

colorbar; 

delete ('DFSM result s.mat'); 


newname = ['DFSM ' filename]; 

cd('C:\Documents and Settings\afbma\Tese\Result Figures'); 
saveas (gcf, newname, 'jpg'); 

cd('C:\Documents and Settings\aflima\Tese\cyclo'); 
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APPENDIX B. FSK/PSK GENERATION CODES 


% FREQUENCY AND PHASE SHIFT COMBINED CODE using Costas Frequency modulation 
% Developed by Antonio Lima and Harsha Tummala 
% July, 2002 

% Phase code for FSK/PSK from IEEE International Radar Conference Paper - 1990 
% “The ambiguity properties of FSK/PSK signals” 

% by J. Patrick Donohoe, and Franklin M. Ingels 


clear all; 
clc; 

disp('FREQUENCY AND PHASE SHIFT CODE (FSK/PSK) USING COSTAS FREQUENCIES******') 


dispC 




') 


% DEFAULT VARIABLES 


A=l; 

fs=15e3; 

SNRdB = 0; 

cpf=l; 

scale=5; 

barker = 5; 

j=sqrt(-l); 


% Amplitude of CW 
% Sample Rate 

% Signal to Noise Ratio 
%Cycles per frequency (> 10) 

% Scaling for plotting time domain graphs 
% Number of Barker Bots for phase modulation 

%j 


% NEW INPUT 
newvar = 1; 
while newvar == 1; 
dispC ') 

disp('WHICH PARAMETER DO YOU WANT TO SET ? ') 
disp('') 

fprintf('l. Amplitude of frequencies - A= %g.\n', A) 
tprintf('2. Sampling frequency - fs (Hz)= %g.\n', fs) 
t^rintf('3. Signal to noise ratio - SNR dB (dB) = %g.\n', SNR dB) 
tprintf('4. Cycles per frequency = %gAn', cpf) 

tprintf('5. Number of bits per Barker code for phase modulation - barker (13/11/7/5)= %g.\n', barker) 
tprintf('6. No changes\n') 
dispC ') 

option= input('Select a option: '); 

switch option 
case 1 

A=input('New amplitude of the carrier signaU '); 
case 2 

fs=input('New sampling frequency (Hz)= '); 
case 3 

SNR_dB=input('New signal to noise ratio (dB)= '); 
case 4 

cpf=input('New number of cycles per frequency = '); 
case 5 

barker=input('New number of bits for Barker Code = '); 
case 6 

newvar = 0; 
end 
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clc; 

end 

% FREQUENCY CHOICES 
newvar = 1; 
while newvar == 1; 
dispC ') 

disp('WHICH COSTAS FREQUENCY SEQUENCE WOULD YOU LIKE TO USE ? ') 
disp('') 

disp(T. 4, 7, 1, 6, 5, 2, 3'); 

disp('2. 2, 6, 3, 8, 7, 5, 1 '); 

disp('3. 4, 1, 6, 7, 5, 3, 2 '); 

disp('4. 5, 6, 8, 2, 4, 3, 1 '); 

disp('5. 4, 9, 2, 5, 1,9,8'); 
disp('6. 2,4,1,1,7,11,9'); 
disp('7. 2, 2, 5, 1, 3, 7, 4 '); 
dispC') 

option2= input('Select an option: '); 

freq=[4 7 1 6 5 2 3; 

2 6 3 8 7 5 1; 

4 1 6 7 5 3 2; 

5 6 8 24 3 1; 

492 5 1 9 8; 

24117119; 

22 5 1 3 7 4]*1000 

switch option2 
case 1 

seq=freq(l,:); 
case 2 

seq=freq(2,:); 
case 3 

seq=freq(3,:); 
case 4 

seq=freq(4,:); 
case 5 

seq=freq(5,:); 
case 6 

seq=freq(6,:); 
case 7 

seq=freq(7,:); 


end 

newvar=0; 

clc; 

end 

% Create spanning Barker frequency 
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minimum=min(seq); 

SAR=ceil(fs/minimum); % Sampling ratio 
tb=l/(fs); % Sampling period 


% This section generates I & Q without COSTAS phase shift and I & Q with Phase shift. The signals are 
generated 

% five times hy the outer loop. The variable 'index' is used to generate a time vector for time domain plots. 
% The signal is generated at seven samples per phase change. 

pw= 1; 

if barker==13 

phase = [ones(l,pw*5),-(ones(l,pw*2)),ones(l,pw*2),-ones(l,pw),ones(l,pw),- 

ones(l,pw),ones(l,pw)];% 13 bits 
elseif barker==l 1 

phase = [ones(l,pw*3),-(ones(l,pw*3)),ones(l,pw),-ones(l,pw),-ones(l,pw),ones(l,pw),-ones(l,pw)];% 
11 bits 
elseif barker==7 

phase = [ones(l,pw*3),-(ones(l,pw*2)),ones(l,pw),-ones(l,pw)];% 7 bits 
elseif barker==5 

phase = [ones(l,pw*3),-(ones(l,pw)),ones(l,pw)];% 5 bits 
end 

index=0; 
for pp = 1:5 

for xx=l:7 

for nn = 1 :barker 

for n=l:SAR*cpf 

IW 0(index+1 )=A*cos(2 *pi * seq(xx) * (n-1) * tb); 

QWO(index+l)=A*sin(2*pi*seq(xx)*(n-l)*tb); 

I(index+1 )=A*cos(2 *pi * seq(xx) * (n-1) * tb) *phase(nn); 

Q(index+1 )=A* sin(2 *pi* seq(xx) * (n-1) * tb) *phase(nn); 
time(index+l)=index*tb; %time vector cumulation 
index = index +1; 
end 
end 
end 
end 

%Power Spectral Density for I with phase shift & with WGN with Signal to noise ratios (SNR) = [0,- 
5,5,10,-10,-20] 
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%for loop makes calculations and plots for each value of SNR for WGN 
[a,b]=size(I); 

SNR=10^(SNR_dB/10); 

power=10*logl0(A^2/(2*SNR));%calculate SNR indB for WGN function 
noise=wgn(a,b,power);%calculate noise at specified SNR 
IN=I+noise; %add noise to I with FSK/PSK phase shift 

IPWON=I; %I with phase shift without noise 

QN=Q+noise; %add noise to Q with FSK/PSK phase shift 

QPWON=Q; %Q with phase shift without noise 

ffs = fs/1000; 

%PLOTS 

dispC') 

pit = input('Do you want to generate plots of the signal (Y/y or N/n) ?Vs'); 
dispC ') 

if(pit=='Y')|(pit==y) 
dispC') 

%Plot Power Spectral Density for I without phase shift 
% figure ; % open new figure for plot 

% psd(IWO,[],fs); %Power Spectral Density of I without Phase shift 
% title(['PSD of I without Phase Shift or Noise']); 


%Plot PSD of 1+ FSK/PSK Phase + WGN and Time Domain of I + FSK/PSK Phase 
figure ;% open new figure for plot 
signal = I+j*Q; 

psd(signal,[],fs);%plot PSD for specified noise SNR 

title(['PSD of FSK-PSK-C-' num2str(option2) num2str(ffs) num2str(barker) num2str(cpf) '-s 

i+j*Q']); 

figure ;% open new figure for plot 
signal = IN+j*QN; 

psd(signal,[],fs);%plot PSD for specified noise SNR 

title(['PSD of FSK-PSK-C-' num2str(option2) '-' num2str(ffs) '-' num2str(barker) '-' num2str(cpf) '-' 
num2str(SNR_dB)' I+j*Q']); 

% %plot time domain signal I with FSK/PSK phase shift and WGN at specified SNR 
% figure ;%open new figure for plot 

% plot(time(l :floor(size(time,2)/scale)),I(l :floor(size(time,2)/scale))); 

% title(['FSK-PSK-C-' num2str(option2)'-' num2str(ffs)'-' num2str(barker) '-s Time Domain']); 

% xlabel('{\itTime - Seconds} '); 

% ylabel('Amplitude'); 

% grid on; 

% 

% Now check to see if signal is correct by plotting phase shift alone and then determining phase shift 
from I+jQ. 

% To determine phase shift, look at the phase angle of I+iO at every 7th time interval. Expect to see the 
FSK/PSK phase 

% function plot repeated 5 times after unwrapping and detrending the phase angle. 

figure;%open new figure for plot 
plot(phase); 
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title(['FSK-PSK-C-' num2str(option2) num2str(ffs) num2str(barker) num2str(cpf) '-s Phase Shift 
for each Costas frequency hop']); 
xlabel('i - index for phase change'); 
ylabel('FSK/PSK Phase Shift - Theta'); 
grid on; 
else 

disp('Signal not plotted') 

^rintf('\n\n') 

end 

% This section generates the fdes for analysis 

INP=IN';%transpose I with noise and FSK/PSK phase shift for text fde 
QNP=QN';%transpose Q with noise and FSK/PSK phase shift for text fde 
IPWONT=IPWON';%transpose I with phase without noise for text file 
QPWONT=QPWON';%transpose Q with phase without noise for text fde 

% % save results in data fdes 

I=INP(:,1); 

Q=QNP(:,1); 

11= IPWONT(:,l); 

QQ=QPWONT(:,l); 

dispC') 

saveresult = input('Do you want to save the new signal (Y/y or N/n) ?','s'); 

if (saveresult == 'Y') | (saveresult =='y') 
ffs=floor(fs/le3); 

save(['FSK_PSK_C_' num2str(option2) '_' num2str(ffs) '_' num2str(barker) '_' num2str(cpf) '_' 
num2str(SNR_dB)],'r,'Q'); 

I=II; 

Q=QQ; 

save(['FSK_PSK_C_' num2str(option2)'_' num2str(ffs)'_' num2str(barker)'_' num2str(cpf) '_s'],'r,'Q'); 
dispC '); 

disp(['Signal and noise save as : FSK PSK C ' num2str(option2)'_' num2str(ffs)'_' num2str(barker)'_' 
num2str(cpf)'_' num2str(SNR_dB)]); 

disp(['Signal only save as : FSK PSK C ' num2str(option2) '_' num2str(ffs) '_' num2str(barker) '_' 

num2str(cpf) '_s']); 

disp(['Directory: ' num2str(cd)]); 

else 
dispC') 

disp('Signal not saved') 

^rintf('\n\n') 

end 
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% FREQUENCY AND PHASE SHIFT CODE (FSK/PSK) for a target frequency distribution 
% Developed by Antonio Lima and Harsha Tummala 
% July, 2002 

% Phase code for FSK/PSK from IEEE International Radar Conference Paper 
% MATCHED FSK/PSK RADAR 

% by B. Jeffrey Skinner, J. Patrick Donohoe, and Franklin M. Ingels 


clear all; 
clc; 

disp( ***********************************************************************^^ 
jjgp(,H=************pREQUENCY AND PHASE SHIFT CODE (PSK/PSK)*************'); 
disp( ***********************************************************************^^ 


%DEFAULT VARIABLES 


A=l; 

fs =15000; 
SNRdB = 0; 
scale=20; 
j=sqrt(-l); 
global N; 
N=128; 
cpp = 5; 


% Amplitude of CW 
% Sampling Frequency 
% Signal to Noise Ratio 
% Scaling for plotting time domain graphs 

%j 

% Number of fsk/psk sections (frequency hops and phase changes) 
% Number of cycles per phase 


% NEW INPUT 
newvar = 1; 
while newvar == 1; 
dispC ') 

disp('WHICH PARAMETER DO YOU WANT TO SET ? ') 
dispC') 

fprintf(T. Amplitude of the carrier signal - A= %g.\n'. A) 
fprintf('2. Sampling frequency - fs (Hz)= %g.\n', fs) 
fprintf('3. Signal to noise ratio - SNR dB (dB) = %g.\n', SNR dB) 
fprintf('4. Number of phase/frequency hops - N = %g.\n', N) 
fprintf('5. Number of cycles per phase - cpp = %g.\n', cpp) 
fprintf('6. No changes\n') 
dispC ') 

option= input('Select a option: '); 


switch option 
case 1 

A=input('New amplitude of the carrier signaU '); 
case 2 

fs=input('New sampling frequency (Hz)= '); 
case 3 

SNR_dB=input('New signal to noise ratio (dB)= '); 
case 4 

N=input('New number of phase/frequency hops ='); 
case 5 

cpp=input('New number of cycles per phase='); 
case 6 

newvar = 0; 
end 
clc; 
end 


airplane; % Compute frequency hops values according to an airplane's frequency response 
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syn test; % Generates frequeney distribution aceording to previous probability distribution 


global target; % Compute Frequency, first call up airplane variable 
global detection; % 

tb=l/(fs); % Sampling period 

SAR=ceil(fs/min(detection)); % Sampling ratio (for smaller frequency) 

phase = rand(l,cpp); %Compute the phase encoding random sequence for each frequency burst 

greater = fmd(phase>=0.5); 

phase(greater) = pi; 

lesser = fmd(phase<0.5); 

phase(lesser) = 0; 


% This section generates I & Q without phase shift and I & Q with Phase shift. The signals are generated 
% five times the number of frequency hops by the outer loop. The variable 'index' is used to generate a 
time vector for time domain plots. 

% The signal is generated at seven samples per phase change. 
index=0; % Time vector for time domain plots. 

for p=l :2% Generate the signal two times and store sequentially in corresponding vectors 
for ii = (1 :N) %Loop for each frequency hop 
f = detection(ii,:); 
for kk = (l:size(target,l)) 
if f == target(kk,l) 

init_phase = target(kk,3); %defmes initial phase of FSK sequence 
end 
end 

for n=l:SAR*cpp %Loop to increment time for single frequency value. 

I(index+l)=A*cos(2*pi*P(n-l)*tb+phase(ceil(n/SAR))+init_phase); %Calculate in phase 
component of signal with phase shift 

IWO(index+l)=A*cos(2*pi*P(n-l)*tb); % Calculate in phase component of signal without phase 
shift 

Q(index+l)=A*sin(2*pi*P(n-l)*tb+phase(ceil(n/SAR))+init_phase); % Calculate quadrature 
component of signal with phase shift 

QWO(index+l)=A*sin(2*pi*P(n-l)*tb); %Calculate quadrature component of signal without 
phase shift 

time(index+l)=index*tb; %time vector cumulation 
index = index +1; 
end 
end 
end 


%Power Spectral Density for I with phase shift & with WGN with Signal to noise ratios (SNR) = [0,- 
5,5,10,-10,-20] 

%for loop makes calculations and plots for each value of SNR for WGN 
[a,b]=size(I); 

SNR=10^(SNR_dB/10); 

power=10*logl0(A^2/(2*SNR));%calculate SNR indB for WGN function 
noise=wgn(a,b,power);%calculate noise at specified SNR 
IN=H-noise; %add noise to I with FSK/PSK phase shift 

IPWON=I; %I with phase shift without noise 

QN=Q-t-noise; %add noise to Q with FSK/PSK phase shift 
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QPWON=Q; %Q with phase shift without noise 

ffs = fs/1000; 

%PLOTS 

dispC') 

pit = input('Do you want to generate plots of the signal (Y/y or N/n) ?Vs'); 
dispC ') 

if(pit=='Y')|(pit==y) 
dispC') 
global range; 
global airplane; 

%Plot Target Frequeney Distribution and Range/Magnitude plot 
figure;% open new figure for plot 
plot(targetC,l), target(:,2)); title('Ship FFT ABS'); grid on 
title(['Original Target Frequeney Probability Distribution']); 
xlabel('Frequeney'); 

ylabel('Normalized Magnitude = Probability'); 

figure;% open new figure for plot 

plot(range, real(airplane(l:64))); title('Ship'); grid on 

title(['Original Target Range / Magnitude Plot']); 

xlabel('Range (ft)'); 

ylabel('Magnitude'); 

% Plot original frequeney distribution histogram and frequeney random firing distribution 

figure;% open new figure for plot 

orient tall; 

subplot(2,l,l), 

hist(deteetion(:, 1 ),N); 

xlabel('Deteetion Index'); 

ylabel('Number of Oeeurenees'); 

fidl=['Target SYNTHETIC']; 

title(fidl); 

subplot(2,l,2), 
bar(target(:,l), target(:,2)); 
xlabel('Deteetion Index') 
ylabel('Probabibty') 
fid2=['Target ORIGINAL']; 
title(fid2); 

%Plot Power Speetral Density for I without phase shift 
% figure ; % open new figure for plot 

% psd(IWO,[],fs); %Power Speetral Density of I without Phase shift 
% title(['PSD of I without Phase Shift or Noise']); 

%time domain plot of in phase signal I with phase shift 
% figure ; %open new figure for plot 

% % plot small portion of time domain signal I so that data will fit meaningfully in figure. 

% %floor(size(time,2)/seale) seleets a small sample of the veetors to plot 
% plot (time(l :floor(size(time,2)/seale)),I(l :floor(size(time,2)/seale))); 

% title(['Time Domain of I with Phase Shift & no Noise']); 
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% xlabel('{\itTime - Seconds} '); 

% ylabel('Amplitude'); 

% grid on; 

% 

%Plot PSD of 1+ FSK/PSK Phase + WGN and Time Domain of I + FSK/PSK Phase 
figure ;% open new figure for plot 
signal = I+j*Q; 

psd(signal,[],fs);%plot PSD for specified noise SNR 

title(['PSD of FSK-PSK-T-' num2str(ffs)num2str(N)num2str(cpp) '-s I+j*Q']); 

%plot time domain signal I with FSK/PSK phase shift and WGN at specified SNR 
figure ;%open new figure for plot 

plot(time(l:floor(size(time,2)/scale)),I(l:floor(size(time,2)/scale))); 
title(['FSK-PSK-T-' num2str(ffs)num2str(N)num2str(cpp) '-s Time Domain']); 
xlabelC {\itTime - Seconds} '); 
ylabel('Amplitude'); 
grid on; 

% Now check to see if signal is correct by plotting phase shift alone and then determining phase shift 
from I+jQ. 

% To determine phase shift, look at the phase angle of I+jQ at every 7th time interval. Expect to see the 
FSK/PSK phase 

% function plot repeated 5 times after unwrapping and detrending the phase angle. 

figure;%open new figure for plot 
plot(phase); 

title(['FSK-PSK-T-' num2str(ffs)num2sfr(N)num2sfr(cpp) '-s Phase Shift']); 
xlabel('i - index for phase change'); 
ylabel('FSK/PSK Phase Shift - Theta'); 
grid on; 
else 

disp('Signal not plotted') 

^rintf('\n\n') 

end 

% This section generates the files for analysis 

INP=IN';%transpose I with noise and FSK/PSK phase shift for text file 
QNP=QN';%transpose Q with noise and FSK/PSK phase shift for text file 
IPWONT=IPWON';%franspose I with phase without noise for text file 
QPWONT=QPWON';%transpose Q with phase without noise for text file 

% % save results in data files 

I = INP(:,1); 

Q = QNP(:,1); 

II = IPWONT(:,l); 

QQ = QPWONT(:,l); 

dispC') 

saveresult = input('Do you want to save the new signal (Y/y or N/n) ?','s'); 

if (saveresult == 'Y') | (saveresult =='y') 
ffs=floor(fs/le3); 
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save(['FSK_PSK_T_' num2str(ffs) num2str(N) num2str(cpp) num2str(SNR_dB)],'r,'Q', 
'detection'); 

I=II; 

Q=QQ; 

save(['FSK_PSK_T_' num2str(ffs)'_' num2str(N)'_' num2str(cpp) ' s'],T/Q','detection'); 
dispC'); 

disp(['Signal and noise save as : FSK PSK T ' num2str(ffs) '_' num2str(N) '_' num2str(cpp) '_' 
num2str(SNR_dB)]); 

disp(['Signal only save as : FSK PSK T ' num2str(ffs)'_' num2str(N)'_' num2str(cpp) '_s']); 

disp(['Directory: ' num2str(cd)]); 

else 
dispC') 

disp('Signal not saved') 

^rintf('\n\n') 

end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% This part may be implemented on a different fde 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function airplane = airplane; 

global target; 
global airplane; 
global range; 

% Target range Profile 

airplane = [complex(0,-0.00000032:0.00000001 :-0.00000016), complex(-0.001,-0.00000015), complex(- 
0.002,-0.00000014), complex(-0.003,-0.00000013), complex(-0.0027,-0.00000012), complex(- 
0.0025,-0.00000011),... 

complex(-0.0022,-0.00000010), complex(-0.0020, -0.00000009), complex(-0.0017, -0.00000008), 
complex(-0.0015,-0.00000007), complex(-0.0023, -0.00000006), complex(-0.003, -0.00000005), 
complex(-0.0015, -0.00000004),... 

complex(0,-0.00000003), complex(0,-0.00000002), complex(-0.0015, -0.00000001), complex(- 
0.003,0),complex(0.028, 0.00000001), complex(-0.0075, 0.00000002), complex(0.014, 

0.00000003), complex(0.013,0.00000004), complex(-0.038,0.00000005),... 

complex(-0.003, 0.00000006), complex(-0.004, 0.00000007), complex(-0.005, 0.00000008), 
complex(-0.015, 0.00000009), complex(0.003, 0.0000001), complex(-0.0015, 

0.00000011:0.00000001:0.00000016), complex(0, 0.00000017:0.00000001:0.00000031), 

complex(0,-0.00000032:0.00000001:0.00000031)]; 

% Target Time pro fde 
airplane t = flipdim(airplane, 2); 

airplane t = real(airplane_t).*(-l) -I- j.*(imag(airplane_t)); 

% FFT of the time record fo the target response 
airplanel = fft(airplane_t); 
airplane_2 = abs(fftshift(airplane_l)); 
phairplane = phase(airplane_l); 

freq_step = 64; % Number of frequency steps to be used 

fss = 12e6; 

fs = fss/2000; %Downsized sampling frequency (from 12MFlz to 6000Ftz) 
frequency = fs/2:fs/2/freq_step:fs - 1; 
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airplanes = airplane_2(64:127); 
total = sum(airplane_3); 
airplane_4 = airplane_3/total; 

target = [(frequency)', (airplane_4)', (pli_airplane(l:freq_step))']; 

range_bin = (fss*2e-9)^(-l); 

range = l:range_bin:frec|_step*range_bin-l; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% This part may be implemented on a different fde 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%% FREQUENCY SYNTHESIS ALGORITHM FOR FSK/PSK TARGET 
%% Frequencies are in the vector "detection" 
function syn test = syn test; 

fid=['Ship']; 
global N; 
global target; 
new_rec= target(:,l:2); 

%NO EDITS BELOW 

% PULL OFF THE DECIMAL DIGITS 

[nn,jo]=size(new_rec); % Determine number 
% of detections or rows 

% Pull off the density 
% decimal digits as strings 
“delta”=sprintf('%6.4f,new_rec(:,2)); 
ss=reshape(“delta”, 6, nn); 
ss=ss'; 

% Turn back into a number 
% with each digit a column 
“delta”s(:, [ 1 ])=str2num(ss(:, 3)); 

“delta”s(:,[2])=str2num(ss(:,4)); 

“delta”s(:,[3])=str2num(ss(:,5)); 

“delta”s(:, [4])=str2num(ss(:, 6)); 

%% kk =4; % NUMBER OF DECIMAL DIGITS 

% Define probabilities 
% pi through p4 

p0=0; 

pl=10^(-l) * sum(“delta”s(:,l)); 
p2=10^(-2) * sum(“delta”s0,2)); 
p3=10^(-3) * sum(“delta”s0,3)); 
p4=10^(-4) * sum(“delta”s0,4)); 

% Determine pi values 
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piO= 0; 
pil= 10*pl; 

pi2= (10*pl)+(10^2 * p2); 

pi3= (10*pl)+(10^2 * p2)+(10^3 * p3); 

pi4= (10*pl)+(10^2 * p2)+(10^3 * p3) + (lOM * p4); 

% Determine test values for 
% the uniform random variable 

ptest0=0; 

ptestl=pl; 

ptest2=ptestl+p2; 

ptest3=ptest2+p3; 

ptest4=ptest3+p4; 

% FILL MEMORY LOCATION 
% SET 1 

ok=0; 

for j=l :nn % RANDOM VARIABLE INDEX 
if“delta”sG,l)~=0 

fori=l:”delta”sG,l) 

memG+ok,l)=j; 

end 

ok=i+ok; 

end 

end 

if“delta”sG,l)~=0 

[mem size xx]=size(mem); 
else 

memsize = 0; 

XX = 0; 

end 


% FILL MEMORY LOCATIONS 
% SET 2 

ok=0; 

for j=l :nn % RANDOM VARIABLE INDEX 
if“delta”sG,2)~=0 

for i=l:”delta”sG,2) 

mem(i+ok+mem_size, 1 )=j; 
end 

ok=i+ok; 

end 

end 

[mem size xx]=size(mem); 

% FILL MEMORY LOCATIONS 
% SET 3 

ok=0; 

for j=l :nn % RANDOM VARIABLE INDEX 
if“delta”sG,3)~=0 

for i=l:”delta”sG,3) 

mem(i+ok+mem_size, 1 )=j; 
end 

ok=i+ok; 
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end 


end 


[mem size xx]=size(mem); 

% FILL MEMORY LOCATIONS 
% SET 4 

ok=0; 

for j=I :nn % RANDOM VARIABLE INDEX 
if “delta”s(j,4)~=0 

for i=I:”delta”s(j,4) 

mem(i+ok+mem_size, I )=j; 
end 

ok=i+ok; 

end 

end 

% _ 

% Now that mem is fdled 
% generate the detections 

% nn number of detections 

for gi=I:N 

% uni is the uniform RV 
% pull the decimal digits off 

uni=rand; 

uni_str=sprintf('%6.4f,uni) ; 
sss=reshape(uni_str, 6, 1); 
sss=sss'; 

dl=str2num(sss(:,3)); 
d2=str2num(sss(: ,4)); 
d3=str2num(sss(:, 5)); 
d4=str2num(sss(:, 6)); 

% Test the RV to find correct 
% index. Then generate detection 

% DETECTION 

global detection; 

if uni >= 0 & uni < ptestl 

rv_index=mem(round(d 1+1), 1); 
detection(gi, [ 1 ])=new_rec(rv_index, (1:1)); 
elseif uni >= ptestl & uni < ptest2 
const=pi l-(100*(pl)); 

rv_index=mem(round(dl*10 + d2 + const +1),1); 
detection(gi, [ 1 ])=new_rec(rv_index, (1:1)); 
elseif uni >= ptest2 & uni < ptest3 
const=pi2-( 1000* (p 1 +p2)); 

rv_index=mem(round(dl*100 + d2*10 + d3 + const +1),1); 
detection(gi,[l])=new_rec(rv_index,(l: 1)); 
elseif uni >= ptest3 & uni < ptest4 

const=pi3-(10000* (p 1 +p2+p3)); 

rv_index=mem(round(dl*1000 + d2*100 + d3*10 +d4 + const +1),1); 
detection(gi, [ 1 ])=new_rec(rv_index, (1:1)); 
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end 


% %_PLOTS 

% % 

% figure 

% 

% orient tall; 

% subplot(2,l,l), 

% hist(detection(:,l),N); 

% xlabel('Detection Index'); 

% ylabel('Number of Occurences'); 

% fidl=[fid,' SYNTHETIC']; 

% title(fidl); 

% 

% subplot(2,I,2), 

% bar(new_rec(:,I), new_rec(:,2)); 

% xIabel('Detection Index') 

% ylabel('Probability') 

% fid2=[fid,' ORIGINAL']; 

% title(fid2); 
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%PAF for FREQUENCY AND PHASE SHIFT CODE (FSK/PSK) 

% Developed by Antonio Lima 
% July, 2002 

% Phase code for FSK/PSK from IEEE International Radar Conference Paper 
% “MATCHED FSK/PSK RADAR” 

% by B. Jeffrey Skinner, J. Patrick Donohoe, and Franklin M. Ingels 

clear all; 
clc; close all; 

load FSK_PSK_T_15_256_10_s; % Loads the I and Q for the P4 signal, 1 cycle N^2 = 128 

cpp = 10; 

NN = 256; 


% load FSK_PSK_C_l_15_5_l_s; % Loads the I and Q for the P4 signal, 1 cycle N^2 = 128 

% seq = 1; 

% barker = 5; 

%cpf= 1; 

SignalI=I(l:size(I,l)/2)'; % Signal repeates after 896 intervals, 128 phase codes*l samples/phase*3 

SAR*5 

SignalQ=Q(l:size(Q,l)/2)'; 

SignalQJ=SignalQ.*j; % Make Q the complex part 
ul=SignalI+SignalQJ; % Create signal and reference signal 
u2=SignalI- SignalQ J; 


tan = 1; 

T=l; 

N=l; 

tx = length(u2); 
tb = T/tx; 
vstep = 500; 
t=[0:tx-l]*tau/tx; 
dt = t(2)-t(l); 


% Adjust number of pulses in the train 
% Number of sambles in the time domain 
% Code(chip) period 

% Steps in the freq shift dimension - (600 FOR COSTAS) 
% Steps in the time dimension, positive 


z=[]; % Empty result matrix 

v=linspace(-50/tau,50/tau,vstep); % Width in the freq shift dimension - (TARGET, 50 FOR COSTAS) 
for m= 1 :vstep 
mult=[]; 

mult = exp(j*2*pi*v(m)*t); 
u2ex= u2.*mult; 

u3 = conj(u2ex); % Transfer in order to use built in function 

c = xcorr(u3,ul).*dt; % Using the built in correlation function 

universaU[]; % Multiplying with the universal function 

dummy = (pi*v(m)*T); 

dummy = dummy + (dummy==0)*eps; 

universal = abs(sin(dummy*N)/(N*sin(dummy))); 

e= (c.* universal); 

z=[z;abs(e)]; 

v(m) = v(m)*tb; % Normalizing the freq shift axis 

end 

t=[fliplr(-t),t(2:tx)]; % Creating a negative time axis 

t=(t/(tb*tx)); % Normalizing the time axis 

z = z/max(max(z)); 


151 




% FOR FSKPSKTARGET, USE THE FOLLOWING 

PLOTS%%%%%%%%%%%%%%%%%%%%%% 

figure 

contour(t,v,z),grid on 
colorbar 

title(['FSK/PSK Target Code, Cpp = num2str(cpp)N hops = ' num2str(NN)]) 
xlabel('Normalized \tau') 
ylabel('Normalized \nu') 

% axis tight 

figure 

plot(t,z((vstep+l)/2,:)), grid on 

title('PSK/FSK Target code, Cut along the 0 doppler axis') 
xlabel('Normahzed \tau') 
ylabel('Magnitude') 
axis tight 

figure 

plot(v,z(:,tx)), grid on 

title('PSK/FSK Target code. Cut along the 0 delay axis') 
xlabel('Normahzed \nu') 
ylabel('Magnitude') 
axis tight 


% FOR FSKPSKCOSTAS, USE THE FOLLOWING 

PLOTS%%%%%%%%%%%%%%%%%%%%%%% 

% figure 

% contour(t,v,z),grid on 
% colorbar 

% title(['PAF for a FSK/PSK Costas Code, Costas sequence = ', num2str(seq) ', Barker Bits = ' 
num2str(barker)', Cycles per Phase = ' num2str(cpf)]) 

% xlabel('Normahzed \tau') 

% ylabel('Normahzed \nu') 

% axis tight 
% 

% figure 

% plot(t,z((vstep+l)/2,:)), grid on 

% title('PSK/FSK Costas code. Cut along the 0 doppler axis') 

% xlabel('Normahzed \tau') 

% ylabel('Magnitude') 

% axis tight 

% 

% figure 

% plot(v,z(:,tx)), grid on 

% title('PSK/FSK Costas code. Cut along the 0 delay axis') 

% xlabel('Normahzed \nu') 

% ylabel('Magnitude') 

% axis tight 

% 

% 

% %Doing it again 
% index = 0; 

% for ii = l:(cpf*375):((size(I,l))/2) 

% for pp = 0:74 
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% Signall(index+1) = I(ii+pp)'; 

% SignalQ(index+l) = Q(ii+pp)'; 

% index = index + 1; 

% end 
% end 
% 

% SignalQJ=SignalQ.*j; % Make Q the complex part 
% ul=SignalI+SignalQJ; % Create signal and reference signal 
% u2=SignalI-SignalQJ; 

% 

% tan = 1; 

%T = 1; 

% N = 1; % Adjust number of pulses in the train 

% tx = length(u2); % Number of sambles in the time domain 
% tb = T/tx; % Code(chip) period 

% vstep = 500; % Steps in the freq shift dimension - (600 FOR COSTAS) 

% t=[0:tx-l]*tau/tx; % Steps in the time dimension, positive 
%dt = t(2)-t(l); 

% z=[]; % Empty result matrix 

% v=linspace(-20/tau,20/tau,vstep); % Width in the freq shift dimension - (TARGET, 50 FOR COSTAS) 
% for m= 1 :vstep 


% 

mult=[]; 


% 

mult = exp(j*2*pi*v(m)*t); 


% 

u2ex= u2.*mult; 


% 

u3 = conj(u2ex); 

% Transfer in order to use built in function 

% 

c = xcorr(u3,ul).*dt; 

% Using the built in correlation function 

% 

universaU[]; 

% Multiplying with the universal function 


% dummy = (pi*v(m)*T); 

% dummy = dummy + (dummy==0)*eps; 

% universal = abs(sin(dummy*N)/(N*sin(dummy))); 

% e= (c.* universal); 

% z=[z;abs(e)]; 

% v(m) = v(m)*tb; % Normalizing the freq shift axis 

% end 

% t=[fliplr(-t),t(2:tx)]; % Creating a negative time axis 

% t=(t/(tb*tx)); % Normalizing the time axis 

% z = z/max(max(z)); 

% 

% %%%%%%%%%%%%%% FOR FSK_PSK_COSTAS, USE THE FOLLOWING 

PLOTS%%%%%%%%%%%%%%%%%%%%%%% 

% figure 

% contour(t,v,z),grid on 
% colorbar 

% title(['PAF for a FSK/PSK Costas Code, Costas sequence = ', num2str(seq) ', Barker Bits = ' 
num2str(barker)', Cycles per Phase = ' num2str(cpf)]) 

% xlabel('Normahzed \tau') 

% ylabel('Normahzed \nu') 

% axis tight 
% 

% figure 

% plot(t,z((vstep+l)/2,:)), grid on 

% title('PSK/FSK Costas code. Cut along the 0 doppler axis') 

% xlabel('Normahzed \tau') 

% ylabel('Magnitude') 

% axis tight 

% 
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% figure 

% plot(v,z(:,tx)), grid on 

% title('PSK/FSK Costas code, Cut along the 0 delay axis') 

% xlabel('Normalized \nu') 

% ylabel('Magnitude') 

% axis tight 

% 

% %%%%%%%%%%% FOR OTHER PLOTS, USE THE 

FIGURES%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% % figure 

% % colormap(jet); caxis=[3 1]; 

% % mesh(t,v,z),grid on 
% % colorbar; 

% % title('PAF for a FSK/PSK Code') 

% % xlabel('Normalized \tau') 

% % ylabel('Normalized \nu') 

% % zlabel('Magnitude') 

% % axis tight 
% % 

% 

% % figure 

% % mesh(t,v,z),grid on 
% % view(0,0) 

% % colorbar 

% % title('Side-view towards the delay axis') 

% % xlabel('Normalized \tau') 

% % ylabel('Normalized \nu') 

% % zlabel('Magnitude') 

% % axis tight 
% % 

% % figure 

% % mesh(t,v,z),grid on 
% % view(90,0) 

% % colorbar 

% % title('Side-view towards the doppler axis') 

% % xlabel('Normalized \tau') 

% % ylabel('Normalized \nu') 

% % zlabel('Magnitude') 

% % axis tight 
% % 


FOLLOWING 
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APPENDIX C. LIST OF LPI RADAR SIGNALS ANALYZED 


SIGNAL NAME 

CARRIER 
FREQUENCY 
OR HOPPING 
SEQUENCE 
(kHz) 

BW 

(Hz) 

NPBB/ 

CPP 

NUMBER OF 
BITS / PHASES 
AND 

MODULATION 

PERIOD 

MODULATION TYPE 

B 1 7 11 1 -6.mat 

1 

1000 

1 

11 

BPSK 

B 1 7 11 1 O.mat 

1 

1000 

1 

11 

BPSK 

B 1 7 11 1 s.mat 

1 

1000 

1 

11 

BPSK 

B 1 7 11 5 -6.mat 

1 

200 

5 

11 

BPSK 

B 1 7 11 5 O.mat 

1 

200 

5 

11 

BPSK 

B 1 7 11 5 s.mat 

1 

200 

5 

11 

BPSK 

B 1 7 7 1 -6.mat 

1 

1000 

1 

7 

BPSK 

B 1 7 7 1 O.mat 

1 

1000 

1 

7 

BPSK 

B 1 7 7 1 s.mat 

1 

1000 

1 

7 

BPSK 

B 1 7 7 5 -6.mat 

1 

200 

5 

7 

BPSK 

B 1 7 7 5 O.mat 

1 

200 

5 

7 

BPSK 

B 1 7 7 5 s.mat 

1 

200 

5 

7 

BPSK 

FR 1 7 16 1 -6.mat 

1 

1000 

1 

16 

FRANK 

FR 1 7 16 1 O.mat 

1 

1000 

1 

16 

FRANK 

FR 1 7 16 1 s.mat 

1 

1000 

1 

16 

FRANK 

FR 1 7 16 5 -6.mat 

1 

200 

5 

16 

FRANK 

FR 1 7 16 5 O.mat 

1 

200 

5 

16 

FRANK 

FR 1 7 16 5 s.mat 

1 

200 

5 

64 

FRANK 

FR 1 7 64 1 -6.mat 

1 

1000 

1 

64 

FRANK 

FR 1 7 64 1 O.mat 

1 

1000 

1 

64 

FRANK 

FR 1 7 64 1 s.mat 

1 

1000 

1 

64 

FRANK 

FR 1 7 64 5 -6.mat 

1 

200 

5 

64 

FRANK 

FR 1 7 64 5 O.mat 

1 

200 

5 

64 

FRANK 

FR 1 7 64 5 s.mat 

1 

200 

5 

64 

FRANK 

F 1 7 250 20 -6.mat 

1 

250 

X 

20 

FMCW 

F 1 7 250 20 O.mat 

1 

250 

X 

20 

FMCW 

F 1 7 250 20 s.mat 

1 

250 

X 

20 

FMCW 

F 1 7 250 30 -6.mat 

1 

250 

X 

30 

FMCW 

F 1 7 250 30 O.mat 

1 

250 

X 

30 

FMCW 

F 1 7 250 30 s.mat 

1 

250 

X 

30 

FMCW 

F 1 7 500 20 -6.mat 

1 

500 

X 

20 

FMCW 

F 1 7 500 20 O.mat 

1 

500 

X 

20 

FMCW 

F 1 7 500 20 s.mat 

1 

500 

X 

20 

FMCW 

F 1 7 500 30 -6.mat 

1 

500 

X 

30 

FMCW 

F 1 7 500 30 O.mat 

1 

500 

X 

30 

FMCW 

F 1 7 500 30 s.mat 

1 

500 

X 

30 

FMCW 


Table 49. Test matrix of LPI radar signals analyzed. 
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PI 1 7 16 1 -6.mat 

1 

1000 

1 

16 

PI 

PI 1 7 16 1 O.mat 

1 

1000 

1 

16 

PI 

PI 1 7 16 1 s.mat 

1 

1000 

1 

16 

PI 

PI 1 7 16 5 -6.mat 

1 

200 

5 

16 

PI 

PI 1 7 16 5 O.mat 

1 

200 

5 

16 

PI 

PI 1 7 16 5 s.mat 

1 

200 

5 

16 

PI 

PI 1 7 64 1 -6.mat 

1 

1000 

1 

64 

PI 

PI 1 7 64 1 O.mat 

1 

1000 

1 

64 

PI 

PI 1 7 64 1 s.mat 

1 

1000 

1 

64 

PI 

PI 1 7 64 5 -6.mat 

1 

200 

5 

64 

PI 

PI 1 7 64 5 O.mat 

1 

200 

5 

64 

PI 

PI 1 7 64 5 s.mat 

1 

200 

5 

64 

PI 

P2 1 7 16 1 -6.mat 

1 

1000 

1 

16 

P2 

P2 1 7 16 1 O.mat 

1 

1000 

1 

16 

P2 

P2 1 7 16 1 s.mat 

1 

1000 

1 

16 

P2 

P2 1 7 16 5 -6.mat 

1 

200 

5 

16 

P2 

P2 1 7 16 5 O.mat 

1 

200 

5 

16 

P2 

P2 1 7 16 5 s.mat 

1 

200 

5 

16 

P2 

P2 1 7 64 1 -6.mat 

1 

1000 

1 

64 

P2 

P2 1 7 64 1 O.mat 

1 

1000 

1 

64 

P2 

P2 1 7 64 1 s.mat 

1 

1000 

1 

64 

P2 

P2 1 7 64 5 -6.mat 

1 

200 

5 

64 

P2 

P2 1 7 64 5 O.mat 

1 

200 

5 

64 

P2 

P2 1 7 64 5 s.mat 

1 

200 

5 

64 

P2 

P3 1 7 16 1 -6.mat 

1 

1000 

1 

16 

P3 

P3 1 7 16 1 O.mat 

1 

1000 

1 

16 

P3 

P3 1 7 16 1 s.mat 

1 

1000 

1 

16 

P3 

P3 1 7 16 5 -6.mat 

1 

200 

5 

16 

P3 

P3 1 7 16 5 O.mat 

1 

200 

5 

16 

P3 

P3 1 7 16 5 s.mat 

1 

200 

5 

16 

P3 

P3 1 7 64 1 -6.mat 

1 

1000 

1 

64 

P3 

P3 1 7 64 1 O.mat 

1 

1000 

1 

64 

P3 

P3 1 7 64 1 s.mat 

1 

1000 

1 

64 

P3 

P3 1 7 64 5 -6.mat 

1 

200 

5 

64 

P3 

P3 1 7 64 5 O.mat 

1 

200 

5 

64 

P3 

P3 1 7 64 5 s.mat 

1 

200 

5 

64 

P3 

P4 1 7 16 1 -6.mat 

1 

1000 

1 

16 

P4 

P4 1 7 16 1 O.mat 

1 

1000 

1 

16 

P4 

P4 1 7 16 1 s.mat 

1 

1000 

1 

16 

P4 

P4 1 7 16 5 -6.mat 

1 

200 

5 

16 

P4 

P4 1 7 16 5 O.mat 

1 

200 

5 

16 

P4 

P4 1 7 16 5 s.mat 

1 

200 

5 

16 

P4 

P4 1 7 64 1 -6.mat 

1 

1000 

1 

64 

P4 

P4 1 7 64 1 O.mat 

1 

1000 

1 

64 

P4 


Table 50. Test matrix of LPI radar signals analyzed. 
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P4 I 7 64 I s.mat 

1 

1000 

P4 I 7 64 5 -6.mat 

1 

200 

P4 I 7 64 5 O.mat 

1 

200 

P4 I 7 64 5 s.mat 

1 

200 

C I 15 10 -6.mat 

4, 7, 1, 6, 5, 2, 3 

- 

C I 15 10 O.mat 

4, 7, 1, 6, 5, 2, 3 

- 

C_l_15_I0_s.mat 

4, 7, 1,6, 5, 2,3 

- 

C I 15 20 -6.mat 

4, 7,1,6, 5,2,3 

- 

C_I_I5_20_0.mat 

4, 7, 1,6, 5,2,3 

- 

C I 15 20 s.mat 

4, 7, 1, 6, 5, 2, 3 

- 

C 2 17 10 -6.mat 

2, 6, 3, 8, 7, 5, 1 

- 

C 2 17 10 O.mat 

2, 6, 3, 8, 7, 5, 1 

- 

C 2 17 10 s.mat 

2, 6, 3, 8, 7, 5, 1 

- 

C 2 17 20 -6.mat 

2, 6, 3, 8, 7, 5, 1 

- 

C 2 17 20 O.mat 

2, 6, 3, 8, 7, 5, 1 

- 

C 2 17 20 s.mat 

2, 6, 3, 8, 7, 5, 1 

- 

FSK PSK C I 15 II I O.mat 

4, 7, 1, 6, 5, 2, 3 

1000 

FSK PSK C I 15 II I s.mat 

4, 7, 1, 6, 5, 2, 3 

1000 

FSK PSK C I 15 II 5 O.mat 

4, 7, 1,6, 5, 2, 3 

200 

FSK PSK C I 15 II 5 s.mat 

4, 7, 1, 6, 5, 2, 3 

200 

FSK PSK C I 15 5 I O.mat 

4, 7, 1,6, 5,2,3 

1000 

FSK PSK C I 15 5 I s.mat 

4, 7, 1,6, 5,2,3 

1000 

FSK PSK C I 15 5 5 O.mat 

4, 7, 1,6, 5,2,3 

200 

FSK PSK Cl 15 5 5 s.mat 

4, 7, 1, 6, 5, 2, 3 

200 

FSK_PSK_T_I5_I28_I0_0.mat 

- 

3600 

FSK PSK T 15 128 10 s.mat 

- 

3600 

FSK PSK T 15 128 5 O.mat 

- 

4200 

FSK PSK T 15 128 5 s.mat 

- 

4200 

FSK PSK T 15 256 10 O.mat 

- 

3600 

FSK PSK T 15 256 10 s.mat 

- 

3600 

FSK_PSK_T_I5_256_5_0.mat 

- 

4200 

FSK PSK T 15 256 5 s.mat 

- 

4200 

T I 7 I s.mat 

1000 


T 12 7 2 s.mat 

1000 and 2000 



Table 51. Test matrix of LPI radar signals ai 
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64 

P4 

64 

P4 

64 

P4 

64 

P4 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 
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